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Abstract: In this paper we solve the equations that describe nucleus- nucleus scattering, in high density 
QCD, in the framework of the BFKL Pomeron Calculus. We found that (i) the contribution of short 
distances to the opacity for nucleus-nucleus scattering dies at high energies, (ii) the opacity tends to unity 
"O at high energy, and (iii) the main contribution that survives comes from soft (long distance) processes for 
^ H large values of the impact parameter. The corrections to the opacity £1 (Y, b) = 1 were calculated and 



it turns out that they have a completely different form, namely ( 1 — f2 — > exp f— Const \JY X) than the 

^ opacity that stems from the Balitsky-Kovchegov equation, which is ( 1 — O — > exp (-ConstY 2 )). We 
^ Jreproduce the formula for the nucleus-nucleus cross section that is commonly used in the description of 
nucleus-nucleus scattering, and there is no reason why it should be correct in the Glauber-Gribov approach. 
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1. Introduction 

High density QCD has been developed for the scattering of the dilute system of partons with the dense 
system of partons [1-8]. The typical example of such processes are deep inelastic scattering at low x in 
which one colorless dipole scatters with the dense target, and/or the hadron-nucleus interaction where 
the dilute system of partons in a hadron interacts with the dense partonic system of the target. The 
main non-linear equations that govern these interactions have been found [5-7] and studied in great detail. 
On the other hand, the scattering of the dense system of partons with the dense system of partons has 
been actively studied [9-16] but with limited success, in spite of the fact that this scattering is closely 
related to nucleus-nucleus scattering. The latter is a well known process that has been studied at RHIC 
experimentally, and in which the key property of high density QCD has been observed. 

In this paper we revisit the nucleus-nucleus interaction in the framework of high density QCD. We 
assume that the dense-dense system interacts at high energy with the effective Lagrangian that can be 
obtained from the BFKL Pomeron Calculus. The final form of the effective action for this type of Pomeron 
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Figure 1: The full set of the diagrams that contribute to the scattering amplitude in the kinematic region Eq. ( jl.l| 
or Eq. (pT20l). 



interaction was formulated in Ref. [9] which we will use in this paper. We will discuss the BFKL Pomeron 
Calculus [17, 18] in the next section, and we will show that this approach provides the set of equations that 
describes nucleus- nucleus interactions, which was originally derived in Ref. [9] (see also Ref. [19] and * ). 

We solve the main equations for the nucleus-nucleus interaction, which is possible thanks to the suc- 
cess in solving the nucleus-nucleus interaction in the framework of the BFKL Pomeron Calculus in zero 
transverse dimensions. This Calculus models the real QCD case, if we neglect the fact that the sizes of the 
interacting dipoles could change during the "one-dipole" to "two-dipoles" decay [4]. The equations that 
have been proven in the kinematic region [22,23] are: 

gS A (b)G 31P e^ Y oc g G 3F A\' Z e^ Y « 1; G^e^ « 1 (1.1) 

where gi is the vertex of the Pomeron- nucleon interaction, G 3 jp is the triple Pomeron vertex, SA(b) is the 
number density of nucleons at fixed impact parameter b and Y = ln(l/x) is the rapidity. In this region the 
nucleus-nucleus scattering amplitude reduces to the sum over the class of net diagrams shown in Fig. |l|. 

2. The BFKL Pomeron Calculus 

2.1 Main ingredients of the BFKL Pomeron Calculus 

The main ingredients of the BFKL Pomeron Calculus are the same as in the Pomeron Calculus in zero 
transverse dimensions (see Fig. [I]), including the Green function of the Pomeron and the vertices of the 

* Unfortunately, the set of equations derived in Ref. [9] has not attracted the attention that it deserves, taking into account 
that it is the first theoretical example of treating the dense-dense system of scattering. We can mention only two papers where 
the numerical solutions to the equations have been discussed and some properties of the general solution have been suggested 
(see Refs. [20,21]). Part of these properties are based on the solution for the BFKL Pomeron Calculus in zero transverse 
dimensions, which does not reflect the analytical solution for this problem given in Ref. [22] 
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Pomeron interaction. In the BFKL Pomeron Calculus only the triple Pomeron vertex and the vertex of 
the Pomeron interaction with the nucleon contribute. The first one can be calculated in perturbative QCD 
which is discussed later on, while the Pomeron - nucleon vertex is a pure non-perturbative input. However, 
in order to make all of the calculations more transparent, we will assume the model where the nucleus is 
a bag of A onia, where each of them consists of a heavy quark and antiquark. The typical distances in 
the onium are small and of the order of I/thq, where tuq is the heavy quark mass. Therefore, in this 
oversimplified model we can perform all estimates in the framework of perturbative QCD. 




Figure 2: The graphical form of the BFKL Pomeron Figure 3: The graphical form of the triple Pomeron 
Green function in the coordinate representation vertex in the coordinate representation. 
(Fig. ||(a)) and the amplitude for onium-onium scat- 
tering (Fig. |(b)). 

One of the main ingredients of this approach is the Green function of the BFKL Pomeron G (x%, X2\x' 1 , x' 2 ) 
(see Fig. |2|). The explicit expression for this Green function can be found in Ref. [18]. The onium-onium 
forward scattering amplitude can be easily written using the Green function in the following form; 



A(Y = lns,t = 0) = a 2 s J d 2 b J d 2 r j d 2 R\^ (r)\ 2 \^ {r)\ 2 G (x' l ,x' 2 ;0\x 1 ,x 2 ;Y) (2.1) 

where as is the QCD coupling , ^ is the wave function of the onium, b is the impact parameter which is 
equal to b = ^ [x\ + x 2 + x[ + x' 2 ) and r = x\ — x 2 and R = x[ — x' 2 are the transverse sizes of the two 
scattering colorless dipoles. The Green function has a simple form in this particular case, namely 

d 2 bG(x 1 ,x 2 ;0\x l ,x 2 ;Y) = J — A (0, v) -j== ( J e^ ¥ (2.2) 

The energy levels w(^) are the BFKL eigenvalues 

u(v) = 2 a s (W) - Re if> Q + iv^j = a s (2^(1) - ^ ( 7 ) - V> (1 - 7)) = a sX (7) (2-3) 

where ip(z) = dlnT(z)/dz and T(z) is the Euler gamma function, 7 = i + iv is the anomalous dimension, 
and as = asNc/tr where -/V c is the number of colors. Finally 

a <°'">=ITtW {2A) 
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We do not need to know the explicit expression for the more general case. Instead of this we will discuss 
some key properties of the more generalized Green function. However before this we give the expression for 
the triple Pomeron interaction [9] which can be written in the coordinate representation for the following 
process. Two reggeized gluons with coordinates x' t and x' 2 at rapidity Y[, which we denote as x^l^!'}, 
decays into two gluon pairs {x'(, x'^\Y 2 } and {x^x^IKj'} due to the Pomeron splitting at rapidity Y. This 
Pomeron splitting is shown pictorially in Fig. ||. 

2 ^T / 2 — 2 — 2 {LuG(x' 1 ,x' 2 ;Y;\x 1 ,X2;Y)) x (2.5) 

c J x l2 x 23 x 13 

xG(xi, x 3 ; Y\x", x 3 ; Y 2 ) G(x 3 ,x 2 ; Y\x 3 , x 2 ; Y 3 ) 



where L\ 2 = x V2 p x p 2 with p = — V and x 12 = {x\ — x 2 ) (2-6) 

For further presentation we require some properties of the BFKL Green function, which are listed 
below [18]: 

1. The general definition of the Green function leads to 



(2.7) 



G- 1 (x 1 ,X2;Y\x' 1 ,x' 2 ;Y') = V \p\ (JL + h\ = (JL + H+\ p\p\ 

Hf( Xl ,x 2 ;Y) = g J d 2 x 3 K{x u x 2 \x 3 ) {f( Xl ,x 2 ;Y) - f(x u x 3 ;Y) - f(x 3 ,x 2 ;Y)} (2.8) 
x 2 

K( Xl ,x 2 \x 3 ) = ( 2 - 9 ) 

x 23 x 13 

2. The initial Green function Go that corresponds to the exchange of two gluons, is equal to 

G (x 1 ,x 2 ;Y\x' 1 ,x , 2 ;Y) = n 2 In f ^ ^ | k | f'f (2.10) 



'*'1,2 / 1', 2 / \ 1,2 x l' 2' 



This form of Go has been discussed in Ref. [18]. 

3. It should be stressed that 

v?v|G (xi,x 2 ;n44;^) = (2-ii) 

= (2 vr) 4 [5^ (xi - xi) 5 (2) (x 2 - x' 2 ) + 6® (a* - x' 2 ) (x 2 - x' x ) ) 

4. At high energy the Green function is also the eigenfunction of the operator L12: 

L 12 G(x 1 ,X2;Y\x' 1 ,x' 2 ;Y') = G(x 1 , x 2 ;Y\x[, x' 2 ;Y') « G(xi, x 2 ; Y\x[, x' 2 ; Y'); (2.12) 

A(U, V) 



-A- 



The last equation holds only approximately in the region where v <C 1, but this is the most interesting 
region which is responsible for the high energy asymptotic behavior of the scattering amplitude. 



2.2 Functional integral formulation of the BFKL Pomeron Calculus. 



The theory with the interaction given by Eq. ( |2.2j ) and Eq. ( p.5| ) can be written through the functional 
integral [9] 

= J D<5>D<$> + e s with S = S + S t + S E (2.13) 

where So describes free Pomerons, Sj corresponds to their mutual interaction while Se relates to the 
interaction with the external sources (target and projectile). From Eq. ( |2.5| ) and Eq. ( |2.13| ) it is clear that 

S = J dYdY'd 2 x 1 d 2 x 2 d 2 x' 1 d 2 x' 2 ^ + (x 1 ,x 2 ;Y)G- 1 (x 1 ,x 2 ;Y\x' 1 ,x' 2 ;Y')^(x' 1 ,x' 2 ;Y') (2.14) 

Si = 2vra| I' dyl I ^Wx3 {(w( , iil2 . n) ^(x 1 ,x 3 -Y')^(x 3 ,x 2 ;Y > ) + h.c} (2.15) 
l\ c J J x l2 x 23 x 13 

For Se we have the local interaction both in terms of rapidity and coordinates, namely 

S E = - J dY'd 2 x 1 d 2 x 2 {<S>(x 1 ,x 2 ;Y')T pr (x u x 2 ;Y) + <fr + (xi,x 2 ; Y') T tar (xi,x 2 ; Y)} (2.16) 

where T pr (rt ar ) stands for the projectile and target, respectively. The form of the functions r depend on 
the non- perturb at ive input in our problem. In our simple model of the nucleus they have the following 
form 



= 5{Y- Y') |tf (x 12 ) \ 2 S Al [B-bj; r tar = 8(Y'-0) |* (x 12 ) | 2 S Ai (b) (2.17) 

where Y = Ins and S Ai (b) is the number of nucleons at a given impact parameter b in the nucleus A^. B is 
the impact parameter between the centers of the two nuclei and b is the position of the interacting nucleon 



(onium) in the target. In Eq. ( 2.17 ) we neglect the impact parameter of the onium - onium interaction in 
comparison with the impact parameters of the nucleons (onia) in the nuclei. Indeed, the impact parameter 
of the onium - onium interaction is of the order of the typical onium (nucleon) size, which is much less than 
the nucleus size. This is a key assumption of the Glauber approach [24] which restricts the value of energy 
until we can trust this approach. The radius of the onium-onium interaction increases with energy, and at 
ultra high energy it will be larger than the nucleus size. A discussion of this increase in the framework of 
QCD can be found in Ref. [25]. This assumption means that we can integrate over all impact parameters 
in the Pomeron interaction, and actually all fields $ {xi,x 2 ) and $ + (xi,x 2 ) can be considered to be the 
fields that depend only on the dipole sizes ($ (#12) and & + (^12))- 
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2.3 The equation for nucleus-nucleus scattering 



Fig. U shows the Glauber re-scatterings that give the largest contribution to the nucleus-nucleus scattering 
amplitude. The onium (nucleon) - onium (nucleon) scattering amplitude is of the order of a 2 s and the 
exchange of one BFKL Pomeron leads to the following contribution to the nucleus-nucleus scattering 
amplitude 



Ajp(Y,B) = a 2 s J d 2 bS Al (B-b} S M {b) J d 2 r d 2 R\^ (r)\ 2 \V (R)\ 2 J d 2 b' G [R; 0|r, b'; Y) (2.18) 
\J d 2 bS Al (B-b) S M (b) j d 2 r\V(r)\ 2 g(r;Y) oc a% A 1 / 3 A l J 3 ^(for A x < A 2 ) Y 



a 2 s 



2 /3 

Notice that the factor A± stems from the b integration. The contribution of the triple Pomeron interaction 



leads to the amplitude shown in Fig. |5| Using Eq. (2.5) and the function Q (r;Y) that has been introduced 



in Eq. ( 2.18| ), this contribution can be rewritten in the form 

" ] r d 2 x\ d 2 %2 d 2 x 3 



2a s *% I d 2 bS Al (B - b) S 2 M (6) / dY' [ d X \ d f \ X3 L 12 Q (x 12 ; Y - Y>) (2.19) 

J ^ ' JO J x \2 x 23 x 13 

(x 23 ; Y' - 0) g (x 23 ; Y> - 0) oc (a 2 A 1 / 3 A^ 3 A 2 ' 3 e"® y) U A ^ e -(o) 



Ajp (Y, B) (a 2 s A 2 



2^1/3^(0)7 



Eq. (|2.19|) shows that if 

a 2 A\ /3 Al /3 e^ Y ~ 1; but a 2 s aI /3 e^ Y < 1 (2.20) 

the diagrams with triple Pomeron interactions are small and only the Glauber - type re-scatterings shown 
in Fig. |], give the largest contribution to the nucleus-nucleus scattering amplitude. However, if 

a|<V(°) y ~ 1 while a|e^°) y < 1 (2.21) 

the "net" diagrams of Fig. |l| will also contribute. The set of "net" diagrams (see Fig. [l]) has two topological 
characteristics: (1) these diagrams are two nuclei irreducible or, in other words, any diagram cannot be 
redrawn as two parts separated by two nuclei states in the s-channel; and (2) the smallness of each triple 
Pomeron vertex (G 3 jp oc a|) is compensated by a large factor of A 1 ^ 3 or aIJ 3 . The smallness of a% e w ^°^ y 
leads to the fact that we can neglect the BFKL Pomeron loops or corrections to Pomeron vertices which 
are of the order of / dY 1 dY 2 G z]P {Y 1 ) G 3 jp(y 2 ) exp (w(0)(Yi - Y 2 )) oc a 2 s exp (w(O)Y). 

It is obvious that in order to sum both the Glauber re-scatterings and the net diagrams, we need to 
search for the nucleus-nucleus scattering amplitude in the form^ 



^The proof of this form can be found in Ref. [9] in the framework of the BFKL Pomeron Calculus and in Ref. [22], where 



an explanation can be found of how Eq. (2.25) derives from the parton cascade in the Pomeron Calculus in zero transverse 
dimensions. 
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A(AA;s,B) 



i[l - exp ( -~n(s,B] 



where 



£l(s,B) = Ajp(Y,B) + ... = net diagrams 
a|< 3 < 3 ^ /3 S -(°)sfa|^ /3 8«°\ a 5 A 1/3 ^(°) 



(2.22) 
(2.23) 



(for At < A 2 ) 



Indeed, the Glauber form of Eq. ( 2.22| ) sums all two nuclei reducible diagrams (see Fig. || and Fig. ||) 
while net diagrams, being two nuclei irreducible, provide the contribution VL/2 (see Fig. ||). The normal- 



ization of Eq. ( 2.22 ) is the following 



utot (AA; s) = 2 j d 2 B ImA (AA; s, B) 



(2.24) 



© © 



Figure 4: The Glauber-type re-scattering for the 
nucleus-nucleus scattering amplitude 




Figure 5: The contribution of the triple Pomeron 
vertex to the Glauber-type formula. 



The equations that sum the net diagrams of Fig. [T] can be derived from [9, 19] from the following 
equations of motion 



8S 



ss 



5$ (x,y) 



8<S>+ (x,y) 



0; 



(2.25) 



Using the action of Eq. ( |2 . 1 3| ) , one can derive from Eq. ( 2.25| ) the system of two equations of motion 
for the theory given by Eq. (|2.13| ): 
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- ay, " = (2-26) 
g J d 2 zK(x,y\z) [mx,z;Y')) + (&(z, y; Y')} - {$(x,y;Y')) - ^as{®{z,y;Y')<S>{z,y;Y'))} 



' " 5 4vra 5 J d 2 z ^ G (x, y; Y'\x' , y'; Y>) K (x', y'\z) (z, y> , Y-Y')} d>(x', z; Y>)) 



(2tt) 4 vr 
d{<f> + {x,y;Y -Y' 



BY' 

2ir 



(2.27) 



J d 2 zK (x, y\z) { (*+(*, z-Y- Y')) + ($+(z, y; Y - Y')) - {<f>+(x, y;Y - Y')) - 
4ira s (<S> + (z, y;Y - Y')<5> + {z, y;Y - Y'))} 



(2vr) 4 vr 



1 " 5 4vra 5 J dPz ^f^ Go (x,y;Y'\x',y';Y') K (x',y'\z) ({L zy ,$ (z,y>;Y')} $+(x> ,z;Y - Y')) 



In Eq. ( |2.26| ) and Eq. (2.27) the procedure of averaging is denned as 



D<S>D<S>+0(x,z,Y)e s ^ + } , . 

(0(x,z;Y')) = J— r , v ' ' 2.28 



Eq. ( |2.26| ) and Eq. ( p. 27 ) are general equations for the theory of interacting BFKL Pomerons which 



sum all possible diagrams. For the net diagrams we have the additional property that 
(*(x, z; Y) <S>(z, y; Y)) = (*(x, z; Y)) y; Y)); ($+(x, z; Y) $> + {z, y; Y)) = ($+(x, z; Y)) y; Y)); 

(<£ + (x,z;Y)$(z,y;Y)) = (& + (x,z;Y)) (&(z,y;Y)) (2.29) 



Indeed, drawing all diagrams in the kinematic region of Eq. ( 2.21 ) for ($(x, z;Y) &(z,y;Y)) one can see 



that we have two sets of diagrams, that are disconnected from each other. Introducing 

[D$D$+$(x,z,7)e s [ $ ' $+ l 
N(x,y;Y) = —Anas — ? ; — „ r , and 

N + (x,y;Y) = -A^a s - — = , V gr ' , 2.30 



and using Eq. (p.29|) we rewrite Eq. Q2.26| ) and Eq. Q2.27D in the form 
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dN(x,y;Y>) 



dY' 



(2.31) 



as 
2ir 



d 2 zK(x,y\z) [N{x,z-Y') + N(z,y;Y') - N(x,y;Y') - N(z,y;Y') N(z,y;Y')} 



1 as 
(2vr) 4 vr 



d 2 x'd 2 



if I d%z < /_ L G Q (x,y-Y'\x',y';Y') K(x',y'\z) {L zy , N + (z,y> ,Y - Y')} N(x> , z;Y') 

\ x y ) 



dN+(x,y;Y -Y') 
dY 1 



^ / d 2 zK(x 



(2.32) 



.,_ ; - - ,~,y\z) {n+(x,z;Y-Y') + N+(z,y;Y-Y') - N+(x,y;Y -Y') - 
N + (z,y;Y - Y') N + (z,y;Y -Y')} 



1 a s 



(2tt) 4 7T 




1 





Figure 6: The pictorial representation of Eq. (2.31) )(see Fig. |6|-A). Fig. ||-B shows the initial condition for this 
equation. The black blob denotes the triple gluon vertex while the wavy lines describe BFKL Pomerons. 



The graphical form of the first equation (Eq. fl2.31| )) is shown in Fig. |6|-A. The equation for N + has 
the same graphical form with the replacement N —> N + and N + — > N. Note, that the initial condition for 
N + has the form of Fig. |6|-B but with the Pomeron attached to the upper nucleus. 

Using the initial condition of Fig. [6|-B one can check that the iterations of Eq. (2.31) reproduce the set 
of net diagrams of Fig. |l|. Comparing Fig. [l] with Fig. ^ one can see that Eq. Q2.31 ) sums the net diagrams, 
proving that the set of equations of Eq. (|2.3l| ) and Eq. (|2.32|) describe the nucleus-nucleus scattering. 
Eq. fl2.29|) leads to the second and the third terms of equation. These terms reflect the key property of the 
net diagrams, namely that three Pomerons in a triple Pomeron vertex generates the two sets of diagrams 
that are disconnected from each other (see Fig. 0-A, where these two sets of diagrams are shown in bold 
and normal wavy lines). In other words, each Pomeron line that connects two sets of diagrams gives the 
contribution of the order of a 2 s e^°) y < 1 (see Fig. 0-B for example). 

By assuming that N + is small, Eq. fl2.31|) reduces to the familiar Balitsky-Kovchegov (BK) equation. 
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A 





Figure 7: The diagrams for N. In Fig. we denote by bold wavy lines and by normal lines, the Pomerons that 
contribute to the different N and N 2 terms of Eq. (2.31) (see Fig. ^|-A). In Fig. ^-B the Pomeron that connects the 
two sets of diagrams is shown in blue. One can see that its contribution is suppressed by a factor of a 2 s s u ^ . 



The typical example of this type of situation is deep inelastic scattering (DIS), but even in this case 
we have to be careful [26]. Indeed, the BK equation sums the class of "fan" diagrams shown in Fig. [|. 
For large photon virtualities Q, it is sufficient to take into account just the scattering of one dipole whose 
transverse size is of the order of 1/Q. However, the typical sizes of the dipoles inside of the diagrams of 
Fig. |8]can reach values, such that it is no longer justified to use this class of "fan" diagrams [1,26]. 

It is interesting to notice that in the B-K limit, N being the Green function of two gluons can be 
viewed as the amplitude of the dipole scattering. In the case of nucleus-nucleus scattering the scattering 
amplitude of a dipole that interacts with two nuclei is a more complicated observable which for small N 
and N + is equal to N + N + . 

It should be stressed that the proof of Eq. (2.31) and Eq. ( |2.32| ), as well as the relation to BFKL 
Pomeron diagrams were discussed in Refs. [9,19-21] as well as in Ref. [22] for the BFKL Pomeron Calculus 
in zero transverse dimensions. In this section we just introduce notations and the main ingredients of this 
approach for completeness of presentation. 



Y 




Y 



= Y 







Figure 8: The "fan" diagrams that are summed by 
the Balitsky-Kovchegov equation. The zig zag lines 
denote the BFKL Pomeron. 








Figure 9: The graphical form of the BFKL contri- 
bution of Eq. (O). 



- 10 - 



3. Solution to the main equations deeply in the saturation region. 



3.1 Two saturation scales. 

In Ref. [26], it is shown that actually we are dealing with two saturation scales even in the dilute-dense 
scattering system. Indeed, let us consider the BFKL Pomeron contribution to the case of DIS. From the 
complete set of eigenfunctions of the BFKL equation, one can conclude that (see Eq. (p.2j)) 

G (r, R; t = 0, Y) = j d 2 r' G (r, r'; t = 0, Y - Y') G (r, R; t = 0, Y' - 0) (3.1) 
1 r e+oc dv' f ie+ °° dv r d 2 r' ( r 2 \ iv ( r' 2 \ iv ' ( ^ 

- 7mL tL £J£ 0) (*) + «W- 0) } 

Integrating over r' leads to 5{y — v') and provides the relation given in Eq. Q3.1| ), which is the general 
property of the Green function. On the other hand, each of the Green functions in Eq. ( |3.1|) has its own 
saturation momentum. It is well known [1,27-29] that the equation for the saturation scale does not depend 
on the non-linear terms and can be derived from the knowledge of only the linear part of the equation with 
the BFKL kernel. Using the general equation for the saturation scale [1,27-29], one finds two saturation 
scales for the two Green functions 

G(r,r';t = 0,Y-Y') -> In (r 2 /Vf) = -^lA (y - y') ; (3.2) 

1 7cr 



G(r',R;t = 0,Y'-Q) -+ In (r' 2 /R 2 ) = (Y - 0) ; (3.3) 

J- Icr 



where Q 2 = 1/r 2 and 7 = \ + iv while 7 cr can be found from the equation 

du (7) I _ u (7 cr ) 

dj 1 - Icr 

resolving these equations we obtain two saturation scales: 

,2 2\ ^ (7c ^ 



(3.4) 



^(Qlr 2 ) = -^Z(y-y); (3.5) 

-L 7cr 



ln(Ql s R 2 ) = ^M(y'-O); (3.6) 

Icr 



Eq. ( |3.5| ) and Eq. (|3.6| ) can be rewritten in the form 



0?,. = Q 2 expf - (Y-Y')) : (3.7) 

J- 7cr 



Ql, s = Ql lS (y = 0) exp (pO^l (Y> - 0) ) ; (3.8) 

V 1 - 7 cr '7 

It is useful to introduce two scaling variables: 
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z=0 4=2ln(r'Q.(A;Y=0)) 



Figure 10: Two saturation scales for DIS. 



z + = ln(r /2 Qy = Xs Y'+i; z = - In (r' 2 Q\ s ) = X s (Y -¥')-£ (3.9) 

where £ = In (r /2 Q\ S (Y' = 0)) and A s = uj (7 cr ) /(l — 7 cr ). The saturation effects start to be essential 
when both z and z + are positive. For negative z and z + we can safely use the linear (BFKL) evolution 
equation (see Fig. |l0| ). For DIS due to the large value of the photon virtuality, we have the region where 



the size of the dipole is smaller than both of the inverse saturation scales [26,30] (see Fig. 10). However, in 
the case of nucleus-nucleus scattering the situation is quite different. In this case in the entire kinematic 
region of low x we cannot neglect the saturation effects, rather we have to solve the non-linear equations 
(see Fig. 11). It is worth mentioning that 



z + z + = In (r 2 Q 2 (Y)) = ( where Q 2 = Q\ s (Y = 0) e XaY (3.10) 

One recognizes that Q s is the saturation momentum in the Balitsky-Kovchegov equation for DIS. Although 
the idea of two saturation scales has been on the market during the past six years starting with the paper 
of Mueller and Shoshi [26], it still needs some demystifying. For this purpose let us consider deep inelastic 
scattering at low x in the double log approximation (DLA) . The main contribution stems from the kinematic 
region where 

1 > xi > ... > Xi > ... > x; Qo < kx t ± < ... k i>± < ... < k n ^ < Q; (3.11) 
The dipole amplitude in the momentum representation in the DLA looks as follows 



N(Q;Q ;Y) = Const^ exp (2^ a s Y In (QVQ 2 ,)) (3.12) 
and this amplitude satisfies the unitarity constraints N < 1 if 

Q 2 > Q 2 S = Qle msY (3.13) 
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Eq. ( |3.lD can be written in the form 

N (Q; Q ; Y) oc J din (ft^/Qg) AT (Q, Y — Y') iV (fc^, Q ; Y') (3.14) 
The amplitude JV (A^j., Q ; Y') < 1 for 

kt± >Q 2 S = Qle 4 ^y' (3.15) 

while the amplitude N (Q, ki,±; Y -Y') < 1 if 

Q 2 > kl x e^ Y -Y>) (316) 

One can see that the value of ki ± for which both amplitudes are in the region where we can apply pertur- 
bative QCD or, in other words, where both amplitudes are small, has to satisfy the following inequality 

Qls = Q 2 e-^(y-y) > k f > qi = Qi^sv (3 . 17) 



It should be stressed that even when Q 2 is large (Q 2 > Q 2 , see Eq. fl3.13|) ) and the solution of Eq. ( |3.12 ) 



turns out to be small enough to satisfy the unitarity constraints, the partons inside the DGLAP cascade 
with transverse momenta given in Eq. ( p. 17 ), violate unitarity. If Q 2 = Q 2 (see Eq. ( |3 . 1 3|) ) both parts of 



the inequality of Eq. ( 3.17] ) are equal and all partons violate unitarity. For Q 2 < Q 2 we have the situation 



shown in Fig. 11. We hope that this equation clarified the appearance of two saturation scales and explains 



Fig. 10 and Fig. 11 



It should be stressed that Eq. Q3.2| ) as well as the entire discussion in this section, including Fig. [10] and 
Fig. |ll|, are related to the partons (dipoles) that give the main contribution to the total cross section (BFKL 
Pomeron). To illustrate this point, it is instructive to consider partons with k^± >> Q and ki t ± >> Qo, 
for which we expect that saturation effects will be small, which would seem to disagree with Eq. fljU7|) 
*. These partons do not contribute to the deep inelastic scattering in the DLA approach, as we have seen 
in Eq. ( 3.11| ), and therefore, we do not consider them. However, it is worthwhile to discuss why and how 



these parton are not essential. For these partons that we are discussing, Eq. ( |3.1| ) looks like Eq. ( [3.14 ), 
however 

Q 2 



N (Q;k 2 ± ,Y -Y') = Const^- exp (2 \j a s (Y - Y') In (kf \j Q 2 )) (3.18) 



Taking Eq. ( 3.18 ) into account, we can rewrite Eq. ( ]3.1| ) for these partons as 



f Q 2 dln(k 2 ± /Q 2 ) I I 

N (Q; Qo; Y)cx J ±- exp (2 yj a s (Y - Y') In [k 2 JQ^ + 2 yj a s Y< In (k 2 JQ 2 ) J (3.19) 

Integration over k{ t ± in Eq. ( 3.1 9| ) leads to k^±_ — > Q, and hence the resulting contribution to N (Q; Qq\ Y) 



is proportional to exp [2Ja s Y' In (Q 2 /Q 2 ] ) ) < exp (2 J a s Y In (Q 2 /Ql) ) (Eq. (fTID). In the general 



"^We thank our referee, who drew our attention to this point. 
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case of Eq. 



the integration measure changes from d r'/r , to a I r' /vr 2 R 2 . Consequently the integra- 



tion over r' leads to 5 (v — v' — i) instead of 5 (y — v'). Thus the v' integration leads to the contribution 
of such partons that forms a negligible part of the total cross section. It should be mentioned that a much 
more detailed analysis can be found in Ref. [26]. The last item that we wish to mention, is that we are 
discussing the total cross section, whereas the contribution to the inclusive cross section has a completely 
different form, with two saturation momenta, which are different. 



3.2 Solution deeply in the saturation region : linearized equations 

One can see that in the kinematic region where (see Fig. 11) z > and z + > 0, then both amplitudes 
N (r, Ro;Y') and N + (R, r; Y — Y') are deeply in the saturation region. We can find the solution to 
Eq. (|2~31~1) and Eq. (|2~3l ) in this region, using experience of solving the Balitsky-Kovchegov equation [31]. In 
this approach N is replaced with N = 1 + AiV (and similarly for N + ). This, together with the observation 
that constant N + in Eq. ( 2.31 ), and constant N in Eq. ( 2.32| ) don't contribute, we see that the asymptotic 
solution is N = 1. Replacing N with 1 + AN in Eq. ( p.31[) and acting with L xy on both sides, and using 
Eq. ( |2~Til ) leads to; 



d(L xy AN(x,y;Y')) 



dY' 



(3.20) 



~2tt 



d 2 zK(x,y\z) (L xy AN(x,y;Y')) - -± j d 2 z K (x, y\z) (L xz AN+(x, z; Y')) 
In Eq. (gjg) we neglect all contributions of the order (AN) 2 , AN AN + and (AA r+ ) 2 . A similar result is 



found by using the analogous above mentioned approach on Eq. ( 2.32 ). 
Y * z + =0 




z=0 



£=2ln(r'Ck (A;Y'=0)) 



Figure 11: Two saturation scales for nucleus-nucleus scattering. 
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Introducing the new functions n = L xy AN and h + = L xy AN + we see that Eq. ( ^3l|) and Eq. fl2~32T) 
reduce to two linear equations 



dn jx,y,Y') 
dY' 
dn + (x, y, Y') 
dY' 



2vr 



d 2 zK(x,y\z) h(x,y,Y') - -± J d 2 z K (x, y\z) n + [z, y, Y') (3.21) 
■g / d 2 zK(x,y\z)n + (x,y,Y') ~ ^ [ d 2 z K (x, y\z) h [z, y, Y') (3.22) 



We first rewrite these equations in a more convenient form: 



dn, (x, y, Y') 
dY' 



dh + (x, y, Y') 
dY' 



as 
' 2vr 
as 
"2vr 

as_ 
' 2vr 

as 
"2vr 



d 2 zK(x,y\z) n(x,y,Y') - ^ J d 2 z K (x, y\z) n + (x, y, Y') (3.23) 
(x, y\z) |2n + (x,z,y') — n + (x,y,y')| 



d 2 zK 



d 2 zK{x,y\z) h + (x,y,Y') - ^ J d 2 z K (x,y\z) h (x,y,Y') (3.24) 
d 2 zK(x,y\z) hn[x,z,Y') - h{x,y,Y')X 



We calculate J d z K (x,y\z) taking into account that the main contribution in the saturation region 
stems from the decay of the large size dipole into one small size dipole and one large size dipole. Indeed, 



d 2 zK (x\,X2\z) 
27rln(x 2 2 /p 2 ) : 



7T 



X 12^ X 13 



/T-»2 I ,yi 2 

x 13 l x 12 



7T 



X 



131 



K i2 <ixf 3 



+ vr 



2 X 



13 



d\x 



2 
12 



.7; 



13 



Fl2 X 13l 



(3.25) 



2vr^ 



It should be noted that the £ which is defined in Eq. ( 3.25| ), is different from the £ that appears in Fig. [l0| 
and Fig. |ll|. At this stage we have introduced the artificial cutoff at small values of the dipole size (p), but 
the physical cutoff is p = l/Q s , as discussed in Ref. [31]. It should be noticed that for the opposite case, 
when a dipole decays to two dipoles of larger sizes, there is no logarithmic contribution, since 



J d 2 z K (xi, X2\z) — > TT J 



VQ2 dx 2 



13 



,'4 
13 



(3.26) 



In the two scale kinematic region, the size of the dipole is restricted as follows (see Fig. [Tl| and Eq. (3.7) 
and Eq. Q); 

I/O?,. > r' 2 > l/Q\ s (3.27) 

Therefore, the natural choice is p = l/Q 2 2s . Using Eq. (gjH) , then Eq. ( gj|) and Eq. simplify to 

the following; 
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dn (x, y, Y') 
dY' 



dn + (x, y, Y') 
dY' 



-a s z + |n (x, y, Y') + h + (x, y, Y') } 



(x,y\z) ^2h + (x,z,Y') - h + (x,y,Y')j 



-a s z + |n (x, y, Y') + n+ (x, y, Y' 



as 
2tt 



d 2 zK(x,y\z) hfi(x,z,Y') - ri(x,y,Y')} 

We use the double Mellin transform to solve Eq. ( p. 28 ) and Eq. ( 3.29j ), namely 
+ y') - / — I — r, i ^'+(1-7)2+ 



n\z 



e+lco doo f e+lo ° dry . 
■ 2^/ • 2^i n(w ' 7)e 



n 



(3.28) 



(3.29) 



(3.30) 
(3.31) 



Substituting Eq. ( 3.30 ) and Eq. (|3.31|) into Eq. (|3.28j ) and Eq. (|3.29|) , adding and subtracting these two 
equations we obtain for S (u), 7^ = n (u>, 7) + n + (w, 7) and A ^w, 7^ = n (w, 7) — n + (w, 7) the following; 



(w + A s (l- 7 )) A(w l7 ) 



2 a s 



d£ (^7) 
c?7 



(w + A s (l -7)) S(w,7) = w(7)A(w,7) 



w(7)S(a;,7) 



(3.32) 
(3.33) 



where w (7) is given by Eq. (|2.3[). Solving Eq. (3.32) we find 



E(w,7) 



where we define 



s » w » p { - /„ ^ - — — } 



2a s u (7O 

S (7) exp { L (7) + F (7) w + K (7) } 



(3.34) 
(3.35) 



Hi) 
F( 7 ) 

A- (7) 



1 n d 7 ' 

2a 5 7 w(7 ; ) 
as Jo ^(7 ) 

2a s Jo u) (7') 



2a S 



e^y'a; (7') 



(3.36) 
(3.37) 
(3.38) 
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and where So (7) has to be found from the matching with the regions with one saturation scale (see 
Fig. 11). Using Eq. ( |3.34j ) we obtain the solution in the form 



€+ico doj f e+ico d~i „ . , 
T~- I S o (7) e 

I e— too ^ l J e—ioo 



#(w,7,y',* + ) 



s(z+ r) 

J e— too " J e— too " 

* (w, 7, Y", z + ) = oo(Y' + F (7)) + (1 - 7) z + + uj 2 L (7) 
where the function K (7) which appears in the exponential function in Eq. (|3.35| ) has been absorbed 



(3.39) 
(3.40) 



into the function So (7) in Eq. ( |3.3£ ). The integral over uj in Eq. ( |3.39| ) can be taken explicitly leading to 
the following expression 

H*\r) =- r oo ^^h)yi^ex P ((i- 7) 



/ e— loo 

re+ioo ^ _ 

/ ^— s o (7) exp (1 - 7) z + - 

/ e _ ioo 2m 



4L(7) 



(Y' + F( 7 )y 

4L( 7 ) 

(y' + 2F( 7 )) 



(3.41) 
(3.42) 



where passing from Eq. ( 3.41 ) to Eq. ( 3.42 ), a factor of (vr/L (7))^ and exp (— F 2 (7) /4L (7)) have 
n absorbed into the function Sq (7). From Eq. ( 3.32 ) and Eq. (|3.35| ) we can deduce A (a;, 7), viz; 



A (a;, 7) 



(a; + A s (l- 7 )) 
07(7) 



S (7) exp{L( 7 )a; 2 + F(7)w + ^(7)} 



(3.43) 



Using the double Mellin transform of Eq. (3.30) we can transform A (a;, 7) to A (z + ,Y ') which is a 
function of z + and Y' as 



A^.n - r~ £ r ^ +x i\-^ M (3.44, 

2m J e _ io0 2m uj (7) 



e— too 



exp { (1 - 7) z + + (Y' + F (7)) w + L (7) w 2 } 



where a factor of exp (^(7)) has been absorbed into the function So (7). The integration over u can 
be solved analytically, leading to 

A(2+ , n . - p^ o(7)l /^_r^ + M^ ,3.45) 

J e — ioo 



2m yu \L(j)\ 2w( 7 )L(7) ui (7) 

(y' ± F( 7 )) 5 
4L( 7 ) 



-{(-^ + - il i^ L } 



e — loo 



2vri uw '\^ 2w (7)1,(7) £j( 7 ) 
4L,-A — } 
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where in passing from Eq. ( 3.45| ) to Eq. ( 3.46 ), a factor of (vr/L ( 7 ))^ 2 and exp (— F 2 (7) /4L (7)) has 
been absorbed into the function So (7). From the definition E (z + ,Y') = n(z + ,Y') + h + (z + ,Y') and 
A (z + ,Y') = h(z+,Y') - h + (z + ,Y'), then using Eq. ( |3.42[ ) and Eq. ( |3.46| ) we can obtain n(z + ,Y') and 
n + (z+,Y') as 



n 



^ S (7) -(1 + (7,F')) exp ((1-7)^ + ^^(^,7) ) (3.47) 

(3.48) 



|I S (7) -(l-e ( 7 ,W)) exp((l- 7 )z+ + y'r(y',7) 



where we define 



where we have absorbed all factors that depend on 7 in the function So (7), which should be found 
from the boundary conditions. For z + = the solution reduces to 



n z 



0, Y> - Yq) = £ |JtSo (7) \ (l + ^ (7, Y ' - Y ) ) exp ( (K' - F ) r (7, F' - Y ) ) (3.51) 



In Eq. ( 3.51 ) we introduced the initial value of rapidity Y = Yq which was assumed to be equal to 
zero above. Let us find the Green function demanding that h {z + = 0, Y' — Yq) = 5 (Y' — Yq). If we find 
such a solution, then integrating it over Yq with an arbitrary function will lead to any boundary condition. 
Actually, as we know at z + = 0, then n = Const (see Ref. [31]) but we will discuss this condition below in 
more detail. Choosing 

So (7) = (3-52) 
a s L (7) 



we can rewrite the pre-exponential factor in Eq. ( 3.51|) in the form 



where h (7) is another arbitrary function of 7 to be determined by boundary conditions. It is clear 
that the dr (7, Y') /d-) term in Eq. (|3.53 ) leads to 5 (Y 1 — Yq) after integration over 7, since we have the 
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freedom to change the integration variable in Eq. ( 3.51 ) to r = — (Y' — Yq + 2F (7))/(4L (7)). The second 
term as well as h (7) gives the function that falls down as e~ Y ' 3 2 at large Y' . Finally, the solution to 
Eq. ( |3.47| ) with the boundary condition n (z + = 0, Y 1 — Yq) = 5 (Y' — Yq) has the form 



n 



(z+,Y'~Y ) = S r +t °° ^{*tM + h(>y)} exp((l- 7 ). + +(y'-yo)r(7,r'-y )) (3.54) 



where So is a constant with respect to Y' and £ and 

%m = i±(m 

UJ 2a 5 L( 7 ) 4d 7 V^(7) 



(3.55) 



The next step is to satisfy to the initial condition for Y' = Yq. Along this line the initial condition is given 
by Eq. ( |2.17 ), namely, § 



n 



(t,r = Y ) = ^x? 2 ln(x? 2 /^) S A (b) = 2a l C / qU ^1 Yq) ^ = <o(b)& (3-56) 



Note that when Y' = Y , then Eq. (|3.50|) implies that dr(i,Y' = Yq) /dry = -d/d-y{F{^) /2L( 7 )}. 
Therefore, we need to find So from the condition 



n 



n 



■e+ioo j 



1 d fF(j) 



2m \2a s L(-f) 2 dj \L(j) 



+ Mt) W 1 " 7 * = Co(&)fc 



&Y' = Yo)=*ol . £ 



1 



F( 7 ) 



- ; + M7) ) (3.57) 



£ _ JOO 2wi V2a 5 L( 7 ) 4a s L 2 (7) w (7) 2a s o; (7) L (7) 
= Co(6)^ 

At 7 —7- 0, -F (7) — > 2\ S L{^) + O (7 3 ) and therefore, the only singularity I/7 2 stems from the 
1/L( 7 ) term in Eq. flp?]), ( assuming that ^(7) has no such singularity). Therefore, we need to find So 
from the condition 



(t,Y' = Yo) = So I' 



t^7 



= (1-7)C 



27ri 2asL (7) 

Assuming that 7 in Eq. (3.58) will be small we see that (see Eq. ( p.36| )) 



7<1 1 



Iq 2asuj{-y') 
Using Eq. ( ^59|) , one can see that Eq. (|3.57j ) is satisfied if; 

S = (o(b)/2a s 



44 



7 



(3.58) 



(3.59) 



(3.60) 



^We hope that the notation fo (b) that we use below, will not be confused with £ in Eq. (3.1C) 
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Therefore, we can consider h (7) = if we are not interested in the correction of the order of l/£ to 
the initial condition of Eq. ( 3.56| ), Finally, the solution which takes into account both initial and boundary 
conditions takes the form 



n(z+,Y'-Yn 



Co (b) 



4a 2 , 

^"5 Je — loo 



xexp((1 _ 7)2+ _ (y >_ yo) y--^w) ) 



CoW 
4a| 



(3.61) 



In order to find n (z, Y — Y'), we notice that Eq. (3.47) at the rapidity Y — Y' is also a solution to the 



equation. Using the fact that from Eq. (3.9) z + = X S Y — z, then Eq. (3.61) can be recast as 



n [z 



■ Y - r) =^IZ£i^m{ l + , ^ Y - Y,) } 

x M p((i- 7 )(A, y -,)- (y -y) (y -^ M) ) 



(3.62) 



However generally speaking, So in Eq. ( 3.62| ) is not a constant, but rather a function of Y . Using Eq. ( |3.60D 
and by choosing this function to take the form 



Co [B-b 

EoOO = ^o(;B-b) exp(-A s Y) = K ' exp(-A 5 y) 



(3.63) 



we obtain the required solution for n + (z, Y — Y'). Finally 



Co (B — b) re+ioo j -1 

n+(z,Y-Y') = V 9 7 / ^ vt-t(i + 0(l,Y-Y')} (3.64) 

( H \ (v v>\ (Y-Y' + 2F( 1 )) ^ 
x exp( - (1-7)2, -(Y-Y) ^ j 

Using the steepest decent method we can estimate the typical value of 7 in the integral of Eq. ( [3.61 ). 
The saddle point equation looks as follows 

d*(Y',z+,>y) , (Y'-Ypf . - / 2FQygp) V 

— ^ — = " z " L - (7SP) " (y " yo) Ui^y) = (3 - 65) 
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The solution to Eq. ( |3.65 ) leads to larger values of 7sp. Indeed, at large 7 we have 



W ( 7 ) ^ K »S -2a 5 M«); ^(7) 



7=m; kCJ>1 1 
> 



4a| In (k) 



; ^(7) ► 



A, 7 2 



After substituting these expressions into Eq. ( |j6g ) we obtain that 



75P 



' (F'- Y ) 2 a| ln(«) 

z+ - a s (y - y ) /2 



and the solution behaves as 



4a| In (k) 



(3.66) 



(3.67) 



72 



(z + , Y' -Y ) oc if (Y' - Y , z + ) exp ( z + - Z) 

where Z = \J^(Y' - Y f a% 1ii(k) (z+ - A, (Y' - Y ) /2) 



and k 



' (^-^o) 2 «| 

z+ - a s (y - y ) /2 



(3.68) 



where if is a function which varies slowly with Y' — Yq and z + . Using the same approach we obtain for 
n + (z, Y — Y') the following expression 



n [z 



,Y-Y') 



oc H~ 




where 


Z + 


and 





z+) 



(3.69) 



4 (Y - Y' f al In(«+) (-z + A s (Y - Y') /2) 



g ~ ^) 2 5| 
z — A s (Y — Y' ) /2 



Replacing ln(ft) by In (k) one can calculate the pre-exponential factors in Eq. ( |3.6£ ) and Eq. ( |3.69D 
using Eq. ( 3.66|) and the formula 8.432(6) of Ref. [32], namely 



n 



h (Z) = 2((b) In (k) e z+ [a s [Y'- Y 
{Z+) = 2C(B-b)ln(k + ) e-*{a s (Y-Y') 



(z+-\ s (Y'-Y )/2) 



-z + X s (Y-Y')/2) 
Z+ 



K x (Z) + K (Z) } (3.70) 
K X (Z+) + K (Z + ) } (3.71) 



3.3 Solution in the region with one saturation scale 



The solution given by Eq. ( 3.6Sj ), is correct in the region with two saturation scales, or in other words, the 
region between the two lines z + = and z = in Fig. |ll]. In the region to the right of the line z = 0, then 
z + becomes negative and n + is small here. Therefore, the equation for h takes the form 



uh (u, 7) 



-a s 



dh (w, 7) 
S7 



(3.72) 
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instead of Eq. ( 3.32| ) and Eq. (|3.33|) . The solution to this equation is very simple 

n(w,7) = n (7) exp | - a s w 7} 



(3.73) 



where we need to find n (w, 7) from the matching of this solution with the solution in the kinematic region 
with two saturation scales at z = (see Eq. (gjgD). Substituting Eq. ( g7T3|) into Eq. (|3~30| ) we obtain 



n\z 



,Y>) 



e-\-ioo j. , re-\-ioo j_, 

^ / ^1 ^ o(7; y) ^'+(1-7)^- as -7 



(3.74) 



Integrating over w we see that asj = Y' and finally we have 

n (z + , y') = n {Y'; Y) e z+ ~ ( y '" y °) z+ '^ 
which coincides with the solution given in Ref. [33]. At z = 0, then z + = £, 

£ = 



(3.75) 



q| ln((Y'-Y )/\/C-A s Y72) 2 

- — £_x s y'/2 ~~ ^0) ' an< ^ keeping only the factors in the exponent, we obtain 



n (Y';Y) 



*V(Y',Y) 



(3.76) 



with 



V [Y',Y) = -C + (Y 1 - Y )(/a s + \ 



4 in (y - y )/VC-A s (y'-y 



C-A s (y-y ) 

Therefore, the answer in the kinematic region with z < 0, is 

n(z + ,y',y|z+ < 0) = il (smooth function of z+ and Y') e v(Y'Y)+z+-z+(Y'~Y )/ 



(Y> - Y o y (3.77) 



as 



(3.78) 



It is easy to see that the solution for n + (z,Y',Y\z < 0) can be obtained from Eq. ( |3.78j ) with the 
replacement: Y 1 —> Y — Y' and z + — > z which yields; 



n 



(z,Y',Y\z < 0) = H (smooth function of z and Y - Y 1 ) e v + (Y',Y\E q . fca$) + » - z(Y-Y')^- g 
with 



a|ln (y - Y')/y/\ 8 (Y -Y')-Q 

r,+ (y',y) = -c + (y - y>) c/a s + a| ^ — - _ Ag(y _ y/) (X - y) 2 (3.80) 

It should be stressed that the variable £(&) that we use in the definition of z + and z is equal to 

£ (ft) = In (V 2 S Al (b) J d 2 b' Q\ s (proton, Y' = 0; b') ) for z+ (3.81) 
t(M - Pj = In (V 2 Sa 2 (b - Pj j d 2 b' Q\ s (proton, Y - Y 1 = 0; 6') ) for z 
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4. The nucleus-nucleus scattering amplitude deeply in the saturation region 



Recalling that n = L xy AN and n+ = L xy AN+ we first need to find AN and AN + . Using Eq. ( |2.12D 
and Eq. ( |3.61| ) we can see that 

AN (z + , Y') = <§ g £2 MM J^i 1 + °W } 
xe Xp (( 1 - 7 ), + -(^-, ) ^±p)) ) 



7 e _ ioo 2vrz 16 7 2 (1 - 7) 2 L{i) I + J 

(y'-y + 2F ( 7 )) 



4« 2 . 

S 1 ■> t—ioo 



x exp (1 - 7 ) z+ - {¥' - Y ) 



4L( 7 ) 



xexp((1 _ 7) ,_ ( ,_ n (I™M) 



Assuming that large values of 7 contribute to the integrals in Eq. ( [4.1|) and Eq. (|4.2| ) , we can rewrite these 
equations in a more economical format, after using the formulae 3.471(9) and 8.432(6) in Ref. [32]. In 
this approach, after integrating over 7 , then Eq. ( ^4.1[) and Eq. ( [4. 2D lead to the following expressions for 
AN and AN + 

AN (z+X) = 2((b) In(«) { ( * + " Xs ^ ' " ^ /2 Y A 4 (S) (4.3) 
+ 2,s(Y'-Y )( Z+ - X ° iY z - Y ° )/2 ) 5 K 5{ Z)} 



AA+ (*, y — Y') = 2((B - b) In e - { ( - + As J + n/2 ) ' A 4 (2+) (4.4) 

In Eq. ( ^4.3[ ) and Eq. ( |4.4| ) we use the same notation convention that was used in Eq. ( |3.T0| ) and 
Eq. flgTH]) . The obvious way to calculate f2 (see Eq. fl2.22| )) is to find N (z + = z,Y' = Y), since 

fl/2 = J d 2 x 12 T pr N (z + = z,Y' = Y) (4.5) 
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Figure 12: The graphical representation of the equation used for calculating the scattering amplitude using the 
i-channel unitarity. 
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Figure 13: The example of the diagram which is described in terms of Eq. (4.6) 



However, in order to use Eq. (|4.5[ ) one needs to know the value of N near to the saturation scale, whereas 
our approach has been developed inside of the two regions of saturation scales. We use the method related 
to i-channel unitarity (see Fig. |l^) which was suggested and adjusted to high density QCD in Refs. [26,34]. 
In this approach 



Afl (C) /2 = J dzAN (Y', z) AN (Y - Y', z) 



(4.6) 



It should be stressed that the "net diagrams" of Fig. [j] possess the following remarkable property, that 
by cutting one BFKL Pomeron line we do not change the integration over the kinematic variables that 
describe other Pomerons. Therefore, these Pomerons are included in the description given by equation 
Eq. ( |4.3D and Eq. (|4.4|) for N and N + . In Fig. 13 we give an example of the diagram which is described in 
terms of Eq. ( |4,6j ) . This figure illustrates the fact that the topology of the net diagrams (see Fig. ||) is very 
essential to the derivation of Eq. ( f4.6[ ) . The Pomeron loops cannot be taken into account using Eq. (|4.6p . 
The net diagrams have another remarkable feature, namely that all paths in the diagram that start from 
one nucleus, also finish at the other nucleus. In other words, there are no loops in the diagram. Each path 
can be cut at some value of rapidity, and contributes to Eq. (|4.6| ). Keeping only the factor in the exponent 
we obtain that 
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AO (C)/2 



dz exp 



Z-Z^ 



(4.7) 



J dz exp ^ 



4 In (y-y)/ /i C -z 



(y - y) 2 



CY. 



In (Y'/y/z) 



(Y 



/\2 



where z = z + — ^A s y'. Integrating over z using the steepest decent method we obtain that the saddle 
point in the z- integration is z$p = C/4, which leads to 



An (C) /2 = H exp 



a s In 



y/2 



y 



C/4 / VC/4 



(4. 



It is interesting to note that the dependence on Y' disappears, and the approach to the unitarity bound 
is much milder than in the case of the solution to the BK equation. Indeed, in the saturation region with 
two saturation scales An oc exp ^ — Const y/Y^j while for the BK equation An oc exp ^ — ConstY 2 ^ . 

The function H absorbs all the pre-exponential factors that depend on Y and impact parameters. The 
nucleus-nucleus scattering amplitude can be written using Eq. ( ^22|) - Eq. (§J|) in the form 



ImA (Y; B) = 1 - exp ( - J d 2 bd 2 b' (s Al (b) Sa 2 (b - tj - Afi) 



(4.9) 



where b is the impact parameter of the nucleon inside of one nucleus and B is the distance between the 
centers of two nucleons. b' in Eq. ( |4,9| ) is the impact parameter of the dipole inside of the nucleon. As 
we have discussed above, in our approach we integrated over this impact parameter. The typical value of 
b' is our new dimensional parameter, that determines the range of energy which we can reach within this 
approach. Generally speaking the average b' increases with energy, but we assume that < b' > <C Ra- 
This inequality determines the range of energies where we can trust our approach. The energy dependence 
of < b > cannot be found in our approach, since in the framework of perturbative QCD, < b >oc s A . 
However we know (see the Froissart theorem in ref. [35]) that < b' > can increase only logarithmically, but 
the parameters of such an increase will depend crucially on the non-perturbative parameter: the mass of 
the lightest hadron. The best estimate is to insert into Eq. ( [4.9| ) the experimental value of the inelastic 
cross section for the proton-proton interaction. In doing so we obtain 



W (Y; B) = 1 - exp - a in (pp) / d 2 bS Al (b) S A2 [B-b 



(4.10) 



It should be stressed that Eq. fl4.10| ) flows naturally from our approach, whereas however it is certainly 
incorrect in the usual Glauber- Gribov [36] approximation (see Refs. [37]). However, we would like to stress 
that the form factors S Ai are different from the usual nucleus form factor Sa (&)• Indeed, Sa (b) gives the 
number of nucleons that have impact parameter b and it is equal to 



S A (b) 



dz pa (z, b) 



(4.11) 
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where p is the density of nucleons in the nucleus. Sa (b) is equal to 

/R N 
dzp A (z,b) (4.12) 
-Rjsr 

and it characterizes the situation when all nucleons with the impact parameter b interact as one. 
5. Conclusions 

In this paper we consider nucleus-nucleus scattering at high energies in the framework of the BFKL Pomeron 
Calculus, given by Eq. ( [2,13 ), which follows from the direct sum over Feynman diagrams. It turns out 



that for the dilute-dense system scattering, this approach gives the same results as other approaches, such 
as the dipole approach and the JIMWLK equation. Therefore, it seems reasonable to discuss nucleus 
scattering using the BFKL Pomeron Calculus, since other approaches have failed to successfully derive 
the set of equations for the dense-dense scattering amplitude. However, we would like to mention that in 
spite of the fact that we believe that the BFKL Pomeron Calculus describes the high energy interaction 
in QCD, we cannot exclude the opposite. In particular, it is possible that the four- Pomeron interaction 
could contribute at high energies. The most important result of this paper is the statement that at high 
energies, the equations in the BFKL Pomeron Calculus do not contradict either the s-channel unitarity, or 
the property of crossing symmetry. The main properties of our solution for the nucleus-nucleus scattering 
amplitude, can be summarized as follows. 

1. The contribution of short distances to the opacity 0, dies at high energies. 

2. The opacity tends to unity at high energies. 

3. The main contribution that survives, originates from soft (long distance) processes for large values 
of the impact parameter. 

The corrections to the opacity = 1 — AN that stem from short distances, have been discussed in 
this paper and it was shown that they behave differently from the corrections to the Balitsky-Kovchegov 
equation. Indeed, it turns out that AN oc exp (^—Const\^Y^j while for the BK equation, AN oc 

exp (—ConstY 2 ) . 



The most salient result of this paper is the formula of Eq. ( 4.10 ) that describes nucleus- nucleus colli- 
sions. This formula is instructive, especially since in the usual Glauber-Gribov approach, there is no reason 
to expect that this formula works [37]. 

All of these above mentioned results are based on the BFKL Pomeron Calculus . It should be stressed 
that this Calculus and the action of Eq. ( 2.13j ) has been proven only at N c S> 1. However, we should 



emphasize that the equivalence of this approach with other approaches (see refs. [10, 13] for example) 
that originate from the JIMWLK equation [7], have not been proven. It should be stressed that the two 
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strategies reflect two different fundamental features of QCD, namely that the BFKL Pomeron Calculus 
satisfies, by construction, the i-channel unitarity whereas JIMWLK-based approaches are precise from 
the point of view of s-channel unitarity. In general Eq. (2.31) and Eq. ( 2.32| ) give the equations for the 
dense-dense system of scattering, in the mean field approximation which replaces the BK equation, while 
the solution that was found in this paper has the same legacy as the solution to the BK equation deeply 
in the saturation region, given in Ref. [31]. 

However, the mean field approximation cannot be trusted at extremely high energies. Strictly speaking 
starting from Y = Y m where a|exp ud(0)l^ = 1, both solutions are not valid and we need to take into 
account all kind of enhanced diagrams (see Fig. ^-b for example). However, to calculate such diagrams we 
need to find the Green function of the Pomeron in the mean field approximation. Since at Y = Y m these 
solutions give the amplitude which is close to the unitarity boundary , the Green function can be small^ 
show that such a scenario looks very plausible. Therefore, the problem of taking into account of Pomeron 
loops could be not important for the understanding of the phenomena of saturation. 
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1. Introduction 

High density QCD has been developed for the scattering of the dilute system of partons with the dense 
system of partons [1-7]. The typical example of such processes are deep inelastic scattering at low x in 
which one colorless dipole scatters with the dense target, and/or the hadron-nucleus interaction where 
the dilute system of partons in a hadron interacts with the dense partonic system of the target. The 
main non-linear equations that govern these interactions have been found [5-7] and studied in great detail. 
On the other hand, the scattering of the dense system of partons with the dense system of partons has 
been actively studied [8-14] but with limited success, in spite of the fact that this scattering is closely 
related to nucleus-nucleus scattering. The latter is a well known process that has been studied at RHIC 
experimentally, and in which the key property of high density QCD has been observed. 

In this paper we revisit the nucleus-nucleus interaction in the framework of high density QCD. We 
assume that the dense-dense system interacts at high energy with the effective Lagrangian that can be 
obtained from the BFKL Pomeron Calculus. The final form of the effective action for this type of Pomeron 



- 1 - 



Figure 1: The full set of the diagrams that contribute to the scattering amplitude in the kinematic region Eq. ( jl.l| 
or Eq. ((Til). 



interaction was formulated in Ref. [8] which we will use in this paper. We will discuss the BFKL Pomeron 
Calculus [17, 18] in the next section, and we will show that this approach provides the set of equations that 
describes nucleus- nucleus interactions, which was originally derived in Ref. [8] (see also Ref. [19] and * ). 

We solve the main equations for the nucleus-nucleus interaction, which is possible thanks to the suc- 
cess in solving the nucleus-nucleus interaction in the framework of the BFKL Pomeron Calculus in zero 
transverse dimensions. This Calculus models the real QCD case, if we neglect the fact that the sizes of the 
interacting dipoles could change during the one dipole to two dipoles decay [4]. The equations have been 
proven in the kinematic region [20,21] are: 

gS A (b)G 31P e AY oc g G^A 1 / 3 e AY « 1; G 2 3P e AY « 1 (1.1) 

where gi is the vertex of the Pomeron- nucleon interaction, G 3 jp is the triple Pomeron vertex, Sa(P) is the 
number density of nucleons at fixed impact parameter b and Y = ln(l/x) is the rapidity. In this region the 
nucleus-nucleus scattering amplitude reduces to the sum over the class of net diagrams shown in Fig. |l|. 

2. The BFKL Pomeron Calculus 

2.1 Main ingredients of the BFKL Pomeron Calculus 

The main ingredients of the BFKL Pomeron Calculus are the same as in the Pomeron Calculus in zero 
transverse dimensions (see Fig. |T|), including the Green function of the Pomeron and the vertices of the 

* Unfortunately, the set of equations derived in Ref. [8] has not attracted a deserved atention, being the first theoretical 
example of treating dense-dense system scattering. We can mention only two papers where the numerical solutions to the 
equations have discussed and some properties of the general solution have been suggested (see Refs. [15, 16]). Part of these 
properties are based on the solution for the BFKL Pomeron Calculus in zero transverse dimensions which does not reflect the 
analitical solution for this problem given in Ref. [21] 
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Pomeron interaction. In the BFKL Pomeron Calculus only the triple Pomeron vertex and the vertex of 
the Pomeron interaction with the nucleon contribute. The first one can be calculated in perturbative QCD 
which is discussed later on, while the Pomeron - nucleon vertex is a pure non-perturbative input. However, 
in order to make all of the calculations more transparent, we will assume the model where the nucleus is 
a bag of A onium, where each of them consists of a heavy quark and antiquark. The typical distances 
in the onium are small and of the order of I/thq, where tuq is the heavy quark mass. Therefore, in this 
oversimplified model we can perform all estimates in the framework of perturbative QCD. 



G = 



a) 



G = 



b) 



Figure 2: The graphical form of the BFKL Pomeron 
Green function in the coordinate representation 
(Fig. ^j-a) and the amplitude for onium-onium scat- 
tering (Fig. g-b). 




Figure 3: The graphical form of the triple Pomeron 
vertex in the coordinate representation. 



One of the main ingredients of this approach is the Green function of the BFKL Pomeron G (x%, x 2 \x' 1 , x 2 ) 
(see Fig. |2|). The explicit expression for this Green function can be found in Ref. [18]. The onium-onium 
forward scattering amplitude can be easily written using the Green function in the following form; 



A(Y = lns,t = 0) = a 2 s J d 2 b J d 2 r j d 2 R\$> (r)\ 2 \^ {r)\ 2 G (x' l ,x' 2 ;0\x 1 ,x 2 ;Y) 



(2.1) 



where as is the QCD coupling , ^ is the wave function of the onium, b is the impact parameter which is 
equal to b = ^ [x\ + x 2 + x[ + x' 2 ) and r = x\ — x 2 and R = x[ — x' 2 are the transverse sizes of the two 
scattering colorless dipoles. The Green function has a simple form in this particular case, namely 



G(r,R;t = 0,Y) = J d 2 bG (x 1 ,x 2 ;0\x 1 ,x 2 ;Y) = J 



r 2 \ %v 



ie+oo fa 1 



(2.2) 



The energy levels oj{v) are the BFKL eigenvalues 



2 a s ( -Reif>[- + w 



a s (2^(1) - V> (7) - 1> (1 - 7)) = asX (7) (2-3) 



where ip(z) = dlnT(z)/dz and T(z) is the Euler gamma function, 7 = i + iv is the anomalous dimension, 
and as = asNc/tr where -/V c is the number of colors. Finally 



A(0,i/) 



[1 + Au 2 ] 2 



(2.4) 
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We do not need to know the explicit expression for the more general case. Instead of this we will discuss 
some key properties of the more generalized Green function. However before this we give the expression for 
the triple Pomeron interaction [8] which can be written in the coordinate representation for the following 
process, namely two reggeized gluons with coordinates x^ and x 2 at rapidity Y[ which we denote as 
{x{, a^l Y"/} decay into two gluon pairs {x'{, x'^\Y 2 } and {x 2 , 2/3IY3} due to the Pomeron splitting at rapidity 
Y. This Pomeron splitting is shown pictorially in Fig. |3[ 

ly c J x 12 x 23 x 13 

G(x 1 ,x 3 ;Y\x'(, x 3 ; Yl) G(x 3 ,x 2 ; Y\x 3 , x 2 ; F 3 ') 

where 

L12 = r YiP\P2 with p 2 = —V 2 and r\ 2 = ( x i — x 2 ) 2 (2.6) 

For further presentation we need some properties of the BFKL Green function [18]: 
1. The general definition of the Green function leads to 

G- 1 (x 1 ,x 2 ;Y\x' 1 ,x > 2 ;Y') =p*j%(JL+H) = f^ + H+jp^; (2.7) 



with 



Hf{ Xl ,x 2 ;Y) = ^ / d 2 x 3 K(x u x 2 \x 3 ) (f(x 1: x 2 ;Y) - f( Xl ,x 3 ;Y) - f(x 3 ,x 2 ;Y)); (2.8) 



K( Xl ,x 2 \x 3 ) = (2.9) 

x 23 x 13 



2. The initial Green function (Go), which corresponds to exchange of two gluons, is equal to 



2 2 \ / 2 2 

1 I . v \ 2 1 / 1,1' 2,2' \ , / x l,l' x 2,2' 



G (;ri,x 2 ;r|xi,4;y) = vr 2 In ^ ? In ^ 2 ^ (2.10) 

\ x l,2' x l',2 / \ Xl ' 2 X l',2' / 

This form of Go has been discussed in Ref. [18]. 

3. It should be stressed that 

V 2 V 2 G (xi,x 2 ;y|x' 1 ,x 2 ;y) = (2 7r) 4 ^ 2 \x x - x[) 5 (2 \x 2 - x' 2 ) + 5^(x x - x' 2 ) 5^(x 2 - x[) ) (2.11) 

4. At high energy the Green function is also the eigenfunction of operator L\ 2 : 

LnGixuxxYtf^Y') = —^-G{x u x 2 ;Y\x' x ,x' 2 ;Y') « G(x\, x 2 \ Y\x[,x 2 ; Y'); (2.12) 

A{U, V) 

The last equation holds only approximately in the region where v <C 1, but this is the most interesting 
region which is responsible for the high energy asymptotic behavior of the scattering amplitude. 
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2.2 Functional integral formulation of the BFKL Pomeron Calculus. 



The theory with the interaction given by Eq. ( |2.2j ) and Eq. ( p.5| ) can be written through the functional 
integral [8] 

= J D<S>D$ + e s with 5 = S + 5/ + 5 B (2.13) 
where So describes free Pomerons, Sj corresponds to their mutual interaction while Se relates to the 



interaction with the external sources (target and projectile). From Eq. (2.5) and Eq. (|2.5| ) it is clear that 



S = J dY dY' d 2 x 1 d 2 x 2 d 2 x[d 2 x' 2 ^ + (x 1 ,x 2 ;Y)G- 1 (x 1 ,x 2 ;Y\x' 1 ,x' 2 ;Y')^(x' 1 ,x' 2 ;Y') (2.14) 
Sj = ^ [ dY' [ ^W*3 .{(^gfo^y/)) .*+( Xl , XsiY ')* + (z a ,z 2 ;Y f ) + h.c.} (2.15) 

l\ c J J X 12 ^23 x 13 

For S^; we have the local interaction both in terms of rapidity and coordinates, namely, 

S E = - J dY'd 2 x 1 d 2 X2Mx 1 ,x 2 ;Y')T pr (x 1 ,x 2 ;Y) + <S> + {x u x 2] Y')T tar {x u x 2 -Y)} (2.16) 

where T pr (rj ar ) stands for the projectile and target, respectively. The form of the functions r depend on 
the non-perturbative input in our problem. In our simple model of the nucleus they have the following 
form; 

Tpr = 5{Y- Y') |tf (x 12 ) | 2 S Al (b - 6) ; r tar = 5 (Y' - 0) |* (x 12 ) | 2 S A2 (?) (2.17) 

where Y = Ins and Sa (b) is the number of nucleons at a given impact parameter b in the nucleus A. B is 
the impact parameter between the centers of the two nuclei and b is the position of the interacting nucleon 
(onium) in the target. In Eq. (2.17) we neglect the impact parameter of the onium - onium interaction in 



comparison with the impact parameters of the nucleons (onia) in the nuclei. Indeed, the impact parameter 
of the onium - onium interaction is of the order of the typical onium (nucleon) size, which is much less than 
the nucleus size. This is a key assumption of the Glauber approach [22] which restricts the value of energy 
until we can trust this approach. The radius of the onium-onium interaction increases with energy, and at 
ultra high energy it will be larger than the nucleus size. A discussion of this increase in the framework of 
QCD can be found in Ref. [23]. This assumption means that we can integrate over all impact parameters 
in the Pomeron interaction, and actually all fields $ {xi,x 2 ) and $ + (xi,x 2 ) can be considered to be the 
fields that depend only on the dipole sizes (<I> (212) and <£ + (^12))- 

2.3 Equation for nucleus-nucleus scattering 

Fig. U shows the Glauber rescatterings that give the largest contribution to the nucleus-nucleus scattering 
amplitude. The onium(nucleon)-onium(nucleon) scattering amplitude is of the order of a 2 s and the exchange 
of one BFKL Pomeron leads to the following contribution to the nucleus-nucleus scattering amplitude; 
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A F (Y,B) = a 2 s J d 2 bS Al (B-bj S A2 (b) J d 2 rd 2 R\V (r) | 2 |tf (R) \ 2 J d 2 b' G (R; 0\r, b'; Y) 

oc a|< 3 ^ /3 e^°) y (2.18) 

The contribution of the triple Pomeron interaction leads to the amplitude (see Fig. ||) which is of the order 
of a% A{ /3 A 2 J 3 e 2uJ ^ Y . Therefore if 



a 2 s A\ /3 Al /3 e^ Y ~ 1; but a z s A{'°e^ u > r < 1 



2 a WMo)y 



(2.19) 



the Glauber - type re-scattering shown in Fig. gives the largest contribution to the nucleus-nucleus 
scattering amplitude. However, if 



a 2 s A{ /3 e^ Y ~ 1 while a 2 s e^ Y < 1 



(2.20) 



the net diagrams of Fig. [T] will also contribute. It is obvious that in order to sum both the Glauber 
rescatterings and the net diagrams, we need to search for the nucleus- nucleus scattering amplitude in the 
form; 



where 



A(AA;s,B) = i(l - exp («, B) 

= net diagrams 



n(s,B) = Ajp(Y,B) + 
The normalization of Eq. ( |2.21| ) is the following 



a tot (AA;s) = 2 J d 2 B ImA (AA; s, B\Eq. (|2T2l|) ) 



(2.21) 
(2.22) 

(2.23) 



A, 



A, 




Figure 4: The Glauber-type re-scattering for the Figure 5: The contribution of the triple Pomeron 



nucleus-nucleus scattering amplitude 



vertex to the Glauber-type formula. 



The equations that sum the net diagrams of Fig. [l] can be derived from [8, 19] under the condition that 

5S „ 5S 



5$ (x,y) 



0: 



(x,y) 



0; 



(2.24) 
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Using the action of Eq. (|2.13|) , one can derive from Eq. Q2.24 ) the system of two equations that describes 
the dense-dense scattering in the semiclassical approximation: 



dN(x,y;Y') 



dY' 

as 
2vr 
1 a s 



(2.25) 



d 2 zK(x,y\z) (N(x,z;Y') + N(z,y;Y') - N(x,y;Y') - N(z,y;Y') N(z,y;Y')) 



G {x,y;Y'\x',y';Y') ^^K{x',y'\z) d 2 z (L, y N+ {z,y> ' ,Y -Y>)) N{x',z;Y>) 



dN + (x, y;Y — Y') 



dY' 



(2.26) 



— / d 2 zK(x,y\z) (N + (x, z;Y — Y') + N + (z,y;Y — Y') - N+(x,y;Y - Y') 
2vr J 

N + (z,y;Y -Y')N + (z,y;Y -Y')) 



1 as 
(2?r) 4 7T 



^ / Go (x, y; Y'\x', y'; Y>) fff* K (x',y'\z) d 2 z (L z , y ,N (z,y';Y>)) N+ (x> , z;Y - Y>) 

[x y ) 



r 

r 





Figure 6: The pictorial representation of Eq. 
lines describe gluons. 



The black blob denotes the triple gluon vertex while the helix 



In these equations we use the notation N and N + to mean the following. $ (x, y;Y) = N (x, y; Y) /nots 
and $ + (x,y;Y) = iV + (x,y;Y) /nets, since with this normalization convention, N has the meaning that 
it gives the dipole scattering amplitude. The graphical form of the hrst equation (Eq. (|2.25| )) is shown in 
Fig. ||. Comparing Fig. [l] with Fig. |6]one can see that Eq. ( 2.25 ) sums the net diagrams, proving that the 
set of equations of Eq. ( [2.25 ) and Eq. (|2.26|) describe the nucleus- nucleus scattering. 

By assuming that N + is small, Eq. ( |2.25| ) reduces to the familiar Balitsky-Kovchegov (BK) equation. 
The typical example of this type of situation is deep inelastic scattering (DIS), but even in this case we 
have to be careful [24]. Indeed, the BK equation sums the class of 'fan' diagrams shown in Fig. |7[ For large 
photon virtualities Q it is sufficient to take into account just the scattering of one dipole whose transverse 
size is of the order of 1/Q. However, the typical sizes of the dipoles inside of the diagrams of Fig. can 
reach values, such that it is no longer justified to use these types of 'fan' diagrams [1,24]. 
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Y 




= Y 







Figure 7: The 'fan' diagrams that are summed by 
the Balitsky-Kovchegov equation. The zig zag lines 
denote the BFKL Pomeron. 








Figure 8: The graphical form of the BFKL contri- 
bution of Eq. (O). 



3. Solution to the main equations deeply in the saturation region. 
3.1 Two saturation scales. 

In Ref. [24] it is shown that actually we are dealing with two saturation scales even in the dilute-dense 
scattering system. Indeed, let us consider the BFKL Pomeron contribution to the case of DIS. To complete 
of the set of eigenfunctions of the BFKL equation, one can conclude that (see Eq. fl2.2Q) 



G (r, R; t = 0, Y) = j d 2 r' G (r, r'; t = 0, Y - Y') G (r, R; t = 0, Y' - 0) (3.1) 

R2 . exp{u(v)(Y-Y') + W (i/)(r , -0)} 



V^R* 



.e+oo dv , f ie+oo ^ f ^ / r « N » 



lit I 2-7T 



Integrating over r' leads to 5{y — u') and provides the relation given in Eq. ( |3.1| ), which is the general 
property of the Green function. On the other hand each of the Green functions in Eq. ( |3.1| ) has its own 
saturation momentum. It is well known [1,25,26] that the equation for the saturation scale does not depend 
on the non-linear terms and can be derived from the knowledge of only the linear part of the equation with 
the BFKL kernel. Using the general equation for the saturation scale [1,25,26], one finds two saturation 
scales for the two Green functions; 



G(r, r';t = 0,Y-Y') 
G (/, R; t = 0, Y' - 0) 



ln(r 2 /rf) = (Y-Y>); 
w (7cp) 



In (r?/R*) 



1 -7c 



(Y-0); 



where Q 2 S = 1/r 2 and 7 = h + iv while 7 cr can be found from the equation 



duj (7) ! 



I7=7cr 



1 " 1c 



(3.2) 

(3.3) 



(3.4) 
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resolving these equations we obtain two saturation scales: 



ln(Q 2 s i? 2 ) 



l-7c 

tO (7cr) 
1 - 7cr 



(y-y'); 
(r - 0) ; 



Eq. (3J5) and Eq. (|3.q ) can be rewritten in the form 

w (7 cr ) 



Q? )S = Q 2 exp(-^ 
QL = Qls (X = 0) exp 



1 - 7cr 



(y - o) 



(3.5) 
(3.6) 

(3.7) 
(3.8) 



z + =0 




z=0 



4=2ln(r'Q.(A;Y'=0)) 



Figure 9: Two saturation scales for DIS. 
It is useful to introduce two scaling variables: 



ln(r /2 Qy = A s y'+£; 



ln(r /2 Q liS ) = X s (Y-Y')-t 



(3.9) 



where £ = ln(r Q 2 2 (Y' = 0) and X s = ^_^ r ■ The saturation effects start to be essential when both z 
and z + are positive. For negative z and z + we can safely use the linear (BFKL) evolution equation (see 
Fig. ||[). For DIS due to the large value of the photon virtuality, we have the region where the size of 
the dipole is smaller than both of the inverse saturation scales [24] (see Fig. ||). However, in the case of 
nucleus-nucleus scattering the situation is quite different. In this case in the entire kinematic region of low 



x we cannot neglect the saturation effects, rather we have to solve the non-linear equations (see Fig. |10|) , 
It is worth mentioning that 

z + z + = ln(r 2 Q s (Y)) = ( where Q 2 = Q\ s (Y = 0) e XaY (3.10) 

One recognizes that Q s is the saturation momentum in the Balitsky-Kovchegov equation for DIS. 
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3.2 Solution deeply in the saturation region : linearized equations 



One can see that in the kinematic region where (see Fig. 10) z > and z + > both amplitudes N (r, Rq; Y') 
and -/V (R, r;Y — Y') are deeply in the saturation region. We can find the solution to Eq. ( |2.25| ) and 
Eq. ( 2.26| ) in this region, using experience of solving the Balitsky-Kovchegov equation [27]. In this approach 
N is replaced with N = 1 + AN (and similarly for N + ). This, together with the observation that constant 
N + in Eq. ( |2.25|) , and constant N in Eq. ( [2.26|) don't contribute, we see that the asymptotic solution is 
N = 1, and for AN we have the following equation; 



d(L x , y AN{x,y;Y f )) 



BY' 



(3.11) 



g / d 2 zK(x,y\z) (L x>y AN(x,y;Y')) - ^ / d 2 z K (x, y\z) (L x , z AN + (x, z; Y')) 



a s 



In Eq. ( p.ll[ ) we neglect all contributions of the order (AiV) 2 , AN AN + and (AA^ + ) 2 . 
Y * z + =0 




z=0 



^=2ln(r'C^(A;Y'=0)) 



Figure 10: Two saturation scales for nucleus-nucleus scattering. 



Introducing the new functions h = L xy AN and n + = L xy AN + we see that Eq. ( [2.25 ) and Eq. ( p.26| ) 
reduce to two linear equations 



dn (x,y,Y') 
dY' 
dn + (x, y, Y') 
dY' 



d 2 zK{x,y\z)h + {z,yX) \ (3.12) 



~ / d 2 zK (x,y\z) h(x,y,Y') - — 

Z7T J TV 

■g f d 2 zK(x,y\z)h + (x,y,Y') ~ ^ / d 2 zK (x, y\z) n (z, y, Y') . (3.13) 



We first rewrite these equation in a more convenient form: 
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0ftfoy.rO = _g [<F zK (x,y\z) n{x,y,Y>) - g / ^ : /v (,-. ,/j -,)»- (,-.,,. V) CM I) 



«5 
"2vr 



<i 2 z if (x, y|z) |2n + (x,z,y') — n + (x,y,y')| 



^ + y0 = -§p / d 2 ^(x,y|z)n + (x,y,y') - g / d 2 z K (x,y\z) h (x, y, Y') (3.15) 



•** / d 2 z 
2vr 



if(x,y|z) |2n(x,z,y') - n(x,y,y')j; 



We calculate J* d z K (x, y\z) taking into account that the main contribution in the saturation region 
stems from the decay of the large size dipole into one small size dipole and one large size dipole. Indeed, 



d 2 zK( Xl ,x 2 \z) - n j ^ ^ 2 , = vr f" ^ + vr P* 



13 I ^ 12 ^131 -'p 2 X 13 -'P 2 l rC 12 x 13l 



= 2Tr\n{x{ 2 /p 2 ) = 2vr£ (3.16) 

It should be noted that the £ which is defined in Eq. ( |3.16| ), is different from the £ that appears in Fig. || 
and Fig. |10|. At the moment we introduced the artificial cutoff at small values of the dipole size (p), but 
the physical cutoff is p = 1/Q S , as discussed in Ref. [27]. It should be noticed that the opposite case, when 
a dipole decays to two dipoles of larger sizes, there is not a logarithmic contribution, since 

r-l/Ql j 2 

d 2 zK(x u x 2 \z) vr / — ^ (3.17) 

Jz 2 ^13 



In the two scale kinematic region, the size of the dipole is restricted as follows (see Fig. [l(] and Eq. (3.7) 
and Eq. (I 



1/QL > r' 2 > l/Q\ a (3.18) 

Therefore, the natural choice is p 2 = 1/1/ Q 2 S . Using Eq. fl3,16| ), then Eq. ( |3,14| ) and Eq. ( 3.15| ) simplify 
to the following; 



dn (x,y,Y') 

dY> 



dh + (x, y, Y') 
dY' 



-a s z + {n (x, y, Y') + n+ (x, y, Y') } (3.19) 



(x,y\z) |2n + (x,z,Y') - h + (x,y,Y')} 



-a s z + {n (x, y, Y') + n+ (x, y, Y') } (3.20) 



(x,y\z) ^2h(x,z,Y r ) - h(x,y,Y')\ 
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We use the double Mellin transform to solve Eq. ( [3.19 )and Eq. ( |3.20 ), namely 



■;?. 



J e — ioo ^ 7r * J e — ioo * m 



n 



Je- too Je — ioo z7r * 



ajy'+(i- 7 ) 2 + 



(3.21) 
(3.22) 



Substituting Eq. (^2lj) and Eq. ( ^22|) into Eq. (|3.19|) and Eq. (|3i2C|) , and adding and subtracting 
these two equations we obtain for S (cu, 7) = n (w, 7) + n + (cj, 7) and A (10, 7) = n (w, 7) — n + (w, 7) the 
following; 



(w + A a (l-7)) A(w, 7 ) 



(<j,7) 

2 a s — w(7)E(w,7) 



(w + A s (l -7)) S(w,7) = a; (7) A (a;, 7); 



(3.23) 
(3.24) 



where w (7) is given by Eq. (|2.3|). Solving Eq. (3.23) we find 



S (w,7) = S (7) exp 



(l, + A s (1-7 / )) 2 + ^ 2 (7') 



2a s w (7') 

S (7) exp { L (7) + F (7) u + K (7) } 



(3.25) 
(3.26) 



where we define; 



i(7) 
F( 7 ) 

K( 7 ) 



1 H dry' 
2a s J w(7 ; ) 

as Jo ^(7 ) 

-— /V^Vt 

2as 7o w (7') 



2a 5 



c^y'w (7') 



(3.27) 
(3.28) 
(3.29) 



and where So (7) has to be found from the matching with the regions with one saturation scale (see 
Fig. 10). Using Eq. ( |3.25| ) we obtain the solution in the form 



* (u, 7, Y', z+) = co{Y' + F (7)) + (1 - 7) z+ + uj 2 L (7) 



(3.30) 
(3.31) 
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where the function K (7) which appears in the exponential function in Eq. ( 3.26j ) has been absorbed 
into the function So (7) in Eq. ( 3,30| ). The integral over ui in Eq. ( |3.30| ) can be taken explicitly leading to 
the following expression 



€+ioc dj ~ / Y' 

— S (7) exp (I-7) z + 



2vri 



4L( 7 ) 



4L( 7 ) 

(y' + 2F( 7 )) 



(3.32) 
(3.33) 



where passing from Eq. (EO^) to Eq. fl3~33p , a factor of (it/L(j)) 1/2 and exp (-F 2 (7) /4L (7)) has 
been absorbed into the function Sq (7). From Eq. (|3.23 ) and Eq. (|3.26| ) we can deduce A (w 7 ); 



(q; + A a (l-7)) 

w(t) 



S (7) exp {l (7)^ + ^(7)^ + ^(7)} 



(3.34) 



Using the double Mellin transform of Eq. (3.21) we can transform A (a;, 7) to A (z + ,Y ') which is a 
function of z + and Y' as; 



A(*+ Y') 



€+tOO 



dui 



d 7 (u + X s (1 - 7)) 



<— too 27ri j e _j 00 27ri 



a; (7) 



So (7) 



(3.35) 



exp { (1 - 7) z + + (y ' + F (7)) w + L (7) 



where a factor of exp (K (7)) has been absorbed into the function So (7). The integration over uj can 
be solved analytically, leading to; 



A(z + ,Y') 



TT-^O IT) 



2ni 



7T 



Y' + F( 7 ) A, (I-7) 
L( 7 ) V 2w( 7 )L( 7 ) a; (7) 



exp 



{ (1-7) 



4L( 7 ) 



c^7 



y' + F( 7 ) + A s (l- 7 ) 



exp 



{ (1-7). 



2vri S ° (7) V 2w( 7 )L( 7 ) ' ^(~) 

Y' (y / + 2F( 7 ))- 1 



4L( 7 ) 



7 



(3.36) 



(3.37) 



where in passing from Eq. flQ6| ) to Eq. JO?D , a factor of (vr/L (7)) 1/2 and exp (-F 2 (7) /4L (7)) has 
been absorbed into the function So (7). From the definition S (z + ,Y') = n(z + ,Y') + h + (z + ,Y') and 
A(z+,y') = h(z + ,Y') - h + (z+,Y'), then using Eq. (|3.33|) and Eq. (^37| ) we can obtain h (z+, Y') and 
n+ (*+ Y') as 
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n 



(*+r) : 

[ Z \Y>) : 

where we define; 



llEoW^l + ^y')) exp((l- 7 )z+ + y'r(y',7) 

e+ioo 



^E (7)-(l 
2vri w; 2 v 



7,y') exp (l- 1 )z+ + Y'r(Y',j) 



(3.38) 
(3.39) 



0(7,*") 

r( 7 ,y') 



y' + F( 7 ) + A s (l- 7 ) 



2w (7) L (7) oj (7) 
(y' + 2F( 7 )) 
4L( 7 ) 



_ d /l^+f(7) 



(3.40) 
(3.41) 



where we have absorbed all factors that depend on 7 in the function Eg (7), which should be found 
from the boundary conditions. For z + = the solution reduces to; 



n 



^ E (7) 5 (l + ^ (7, y ' - Y ) ) exp ( (F ' - Y ) r (7, Y ' - Y ) (3)42) 



In Eq. (3.42) we introduced the initial value of rapidity Y = Yq which was assumed to be equal to 
zero above. Let us find the Green function demanding that n (z + = 0, Y' — Yq) = 5 (Y' — Yq). If we find 
such a solution, then integrating it over Yq with an arbitrary function will lead to any boundary condition. 
Actually, as we know at z + n = Const (see Ref. [27]) but we will discuss this condition below in more 
detail. Choosing 



So (7) 



a s L (7) 



(3.43) 



we can rewrite the pre-exponential factor in Eq. (3.42) in the form 



So (7)2 (1 + 0(7, Y') 



Er 



1 



+ 



ldT( 7 ,y') 1 d /F( 7 ) 



2a s L (7) 2 d 7 



4d 7 \L(j) 



+ M7) (3-44) 



where h ( 7 ) is another arbitrary function of 7 to be determined by boundary conditions. It is clear 
that the dr ( 7 , Y') /d 7 term in Eq. (|3.44 ) leads to 5 (Y 1 — Yq) after integration over 7, since we have the 
freedom to change the integration variable in Eq. (]3.42j ) to r = — Y ~^ jj^^ ■ The second term as well as 

h (7) gives the function that falls down as e~ Y ' 3 2 at large Y' . Finally, the solution to Eq. ( |3.38| ) with the 
boundary condition n (z + = 0, Y' — Yq) = 5 (Y' — Yq) has the form 
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n\z ' 



Y'-Yo) = S / " "' £ {*tM + hfr)} exp((l- 7 )z++(y'-yo)T(7,r'-y )) (3.45) 



, e _ loo 2iri I d-y 
where So is a constant with respect to Y' and £ and 

{1) 2a s L( 7 ) 4d 7 V^(7) 



+ Mt) 



(3.46) 



The next step is to satisfy to the initial condition for Y' = Yq. Along this line the initial condition is given 
by Eq. ( |2.17 ), namely, t 



n 



&Y' = Y ) = ^x? 2 ln(xL/r 2 ) S A (b) = ^ Q , £® = ^ ^ = C o(b) ^ (3.47) 



Note that when Y ' = Y then from Eq. (jjUlD dr (7, y ' = Y ) /dry = -d/drf {F (7) /2L (7)}. Therefore, 
we need to find Sq from the condition 



n 



re+ioo j 
«/ e—ioo 



1 



1 d /F( 7 ) 



2vrf V 2 as^(7) 2^7 VMl) 
1 F( 7 ) 



+ /i( 7 )) e (1 - 7)e = CoW^ 



(£,y = y ) = s o^_ iQo ^ i2 ( Kir.) iosL^-um ! LVk.NLI-) 

= Co(6)^ 



+ ^ As( '.:7L +M7)) e (i -^ 

(3.48) 



At 7 — > i 7 (7) — > 2\ s L( r y) + O (7 3 ) and, therefore, the only singularity I/7 2 stems from the 
1/L(j) term in Eq. ( |3.48j ), (assuming that ^(7) has no such singularity). Therefore, we need to find So 
from the condition 



c?7 



1 



=,(1-7)5 



2m 2asL (7) 

Assuming that 7 in Eq. (3.49) will be small we see that (see Eq. ( p. 27] )) 

1 



Hi) 



dj - 



/ 2a s w(Y) 
Using Eq. fl3.50|) , one can see that Eq. ( 3.48j ) is satisfied if; 

S = (o(b)/2a s 



7 <1 1 2 

> T-S-7 



4a| 



(3.49) 



(3.50) 



(3.51) 



Therefore, we can consider h (7) = if we are not interested in the correction of the order of l/£ to 
the initial condition of Eq. ( 3,47p . Finally, the solution which takes into account both initial and boundary 
conditions takes the form 



'We hope that the notation £0 (6) that we use below, will not be confused with £ in Eq. (3.1C) 
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11 u 



,Y'-Y ) 



Co (b) 



dj 1 



4a| J t _ ioo 2vri L (7) 



x exp (1 - 7) z+ - (Y' - y c 



{l + e( 7 ,Y')} 

(Y'-Y + 2F (7)) 



(3.52) 



4L( 7 ) 



Co(b) 
4a| 



7 e _ ioo 2vri L (7) I W ' ' / 



In order to find h (z, Y — Y'), we notice that Eq. ( 3.38 ) at the rapidity Y — Y' is also a solution to the 
equation. Using the fact that from Eq. (3J)) z + = \ S Y — z, then Eq. ( 3.52| ) can be recast as; 



n [z 



X exp((l- 7 )(A s F- z) -(F-r)(- y - i " + 2F W> 



(3.53) 



4L( 7 ) 

However generally speaking, So in Eq. ( [3.53 ) is not a constant, but rather a function of Y. By choosing 
this function to take the form 



EoOO = X [Eq. (plD;£-&j exp(-A s y) = 
we obtain the required solution for ra + (z, K — Y'). Finally, 

Co (B - b) ft+ioo ^ 1 



Co [B - b 

2^7 



exp (-X S Y) 



n (z 



,Y-Y>) 



4a| 



2iri L(j) 



{1 + e( 7 ,y-F') } 



(3.54) 



(3.55) 



Using the steepest decent method we can estimate the typical value of 7 in the integral of Eq. fl3.52|) . 
The saddle point equation looks as follows 

7) _ + (y'-lo) 2 r , 



^7 



4L2 (7 



£ 7 (isp) - (Y 1 - Y ) 



2F 



(tsp) V 

4L( 7SP ) y 



The solution to Eq. ( |3.56 ) leads to larger values of 7s , p. Indeed, at large 7 we have 



w( 7 ) ► -2a s ln(/c); L( 7 ) ► 



4a| In (k) 

After substituting these expression into Eq. ( ggg ) we obtain that 



7 „ / s 7=«; «>i A s 7 Z 

; ^(7) > 



4a| ln(«) 



7sp 



' (y-y ) 2 a| ln(«) 

z+ - X s {Y> - Y ) 12 



(3.56) 



(3.57) 



(3.58) 
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and the solution behaves as 

-,Y'-Y ) oc H (Y' - Y , z + ) exp ( z + - Z) (3.59) 
where Z = y/^(Y'- Y ) 2 a| ln(k) (z+ - X s (Y f - Y Q ) /2) 



n z 



and 



(Y'-Y ) 2 a? s 

- \ s (y - Y ) /2 



where if is a function which varies slowly with y' — Yq and z + . Using the same approach we obtain for 
h + (z, Y — Y') the following expression; 



n + (z,Y-Y') oc H + (Y-Y',z) exp ( - z - Z + ) 

where = J A (Y — Y'f a| 1ii(k+) (-z + A s {Y - Y') /2) 



(3.60) 



- A s (y - Y>) /2 



and 



Replacing ln(/c) by In (re) one can calculate the pre-exponential factors in Eq. ( 3.5S ) and Eq. ( |3.60D 
using Eq. ( 3.57|) and the formula 8.432(6) of Ref. [33], namely 



n 



(Z) = 2C(6)ln(«) e* + {a s (Y'-Y ) ( * + Yq) /2) K,{Z) + K (Z)} 



(3.61) 



n+(Z+) = 2((B-b)ln(K+) e- z {a s (Y-Y>) ( - + K J + n ^ Ifr (Z+) + K (2+) } (3.62) 



3.3 Solution in the region with one saturation scale 



The solution given by Eq. ( 3.59| ), is correct in the region with two saturation scales, or in other words, the 
region between two lines z + = and z = in Fig. [[(]. In the region to the right of the line z = 0, then z 
becomes negative and n + is small here. Therefore, the equation for h takes the form; 



uh (u, 7) 



-a s 



dn (u, 7) 
O7 



instead of Eq. (|3.23j ) and Eq. (3.24). The solution to this equation is very simple 

n(w,7) = n (7) exp j - a s to 7} 



(3.63) 



(3.64) 



where we need to find h (u, 7) from the matching of this solution with the solution in the kinematic region 
with two saturation scales at z = (see Eq. fl3~59D). Substituting Eq. ([T64]) into Eq. ( ggjj ) we obtain 



h(z + ,Y') 



^ / n ( 7 ;y) e <^'+(i- 7 )*+~^7 



(3.65) 
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Integrating over to we see that as'y = Y' and finally we have 



n 



{z + ,Y') = n {Y';Y) e z+ ~ (Y '- Y °1 z+ / as 



(3.66) 



which coincides with solution given in Ref. [31]. At z = then z H 
and keeping only factors in the exponent, we obtain 



i% \n((Y'-Yo)/y/<;-\ a Y'/2) 
t-\ s Y'/2 



ho(Y';Y) = e 



V(Y',Y) 



(3.67) 



with 



r,(Y',Y) = -C + (V - r )C/5 s + 



ax 



In (Y> - Y )/y/C-\s(Y>-Y 



C-Xs(Y'-Y 



(Y> - Y o y (3.68) 



Therefore, the answer in the kinematic region with z < 0, is 



n [z 



,Y',Y\z + < 0) = H (smooth function of z+ and Y 1 ) e ^ Y '' Y ^ Eq - <M» 



+ z+ z+(Y'-Y )/a s 



(3.69) 



It is easy to see that the solution for n + (z,Y',Y\z < 0) can be obtained from Eq. ( |3.69| ) with the 
replacement: Y' — > Y — Y' and z + — > z which yields; 



n 



(z,Y',Y\z < 0) = (smooth function of z and Y — Y') e v + (Y',Y\E q . (^n|)) + z - z (Y-Y>)/- s ^ ^ 



with 



,,+ (!", Y) = -C + (y- 1") a as + 



i 



a. 



In (Y - Y')/y/X 8 (Y -Y')-C 



C - \ S (Y-Y>) 



(Y - Y'Y (3.71) 



It should be stressed that the variable £(&) that we use in definition of z + and z is equal to 



£ (&J = In (r /2 S Al (&J y d 2 b' Q\ s (proton, F' = 0; b') J for 
£ (fl - = In (V 2 Sa 2 (-B - t) J d 2 b' Q 2 2 s (proton, Y - Y' = 0; 6') ) for 



(3.72) 



4. The nucleus-nucleus scattering amplitude deeply in the saturation region 



Recalling that h = L xy AN and n+ = L xy AN + we first need to find AN and AiV+. Using Eq. fl2.12|) 
and Eq. ( |3.52| ) we can see that 



AN(z + ,Y'^ 



(4.1) 



((b) 



e+ioo 



4(3!g J £_joo 



x exp (1 - 7 ) z+ - {¥' - Y ) 



dj 1 

2ni I67 2 (1 - 7) 2 L (7) 

(y'-y + 2F ( 7 )) 



hi 1 + e ^} 



4L (7) 



AiV+ (z,Y -Y') 



( [B — b) i-e+ioo 



4a| 



1 

2^1 I67 2 (1 — 7) 2 L (7) 



+ flfry-rO} (4-2) 



Assuming that large values of 7 contribute to the integrals in Eq. and Eq. (|4.2| ), we can rewrite these 
equation in a more economical format, after using the formulae 3.471(9) and 8.432(6) in Ref. [33]. In 



this approach, after integrating over 7 Eq. (4.1) and Eq. ( (4.2|) lead to the following expressions for AN 
and AN + \ 



AN(z + ,Y') = 2C(6)ln(«) e z+ { 



A s (Y'-Y )/2 



+ 2, s (Y'-Y 0)( Z+ -^ Y z - Y ^ 2 )\ 5i Z)} 



(4.3) 



AiV+ (z, Y — Y') = 2((B - b) In (k + ) e~ z { 



-z + X s (Y-Y')/2 



Z+ 



Ki (Z + ) (4.4) 



+ 2as( y-_ n( _£±^zIV2)\, 



(2 + )} 



In Eq. ( [4. 3D and Eq. (^J) we use the same notation convention that was used in Eq. ( |3.6l| ) and 
Eq. (pg) . 
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The obvious way to calculate n (see Eq. ( 2.21[ )) is to find N (z + = z,Y' = Y), since 



n/2 = J d 2 x 12 ,T pr N (z + = z,Y' = Y) (4.5) 



However, in order to use Eq. ( ^4.5[ ) one needs to know the value of N near to the saturation scale, whereas 
our approach has been developed inside of the two scales saturation regions. We use the method related to 
i-channel unitarity (see Fig. |ll]) which was suggested and adjusted to high density QCD in Refs. [24,32]. 
In this approach; 

AO (C)/2 = J dzAN (Y',z) AN (Y -Y',z) (4.6) 

It should be stressed that the 'net diagrams' of Fig. |l] possess the following remarkable property, that 
by cutting one BFKL Pomeron line we do not change the integration over the kinematic variables that 
describe other Pomerons. Therefore, these Pomerons are included in the description given by equation 
Eq. ( |4.3| ) and Eq. Q4.4j ) for N and N + . In Fig. ^ we give an example of the diagram which is described 
in terms of Eq. (4.6). This figure illustrates the fact that the topology of net diagrams (see Fig. Ill) is very 



essential in the derivation of Eq. ( |4.6| ) . The Pomeron loops cannot be taken into account using Eq. 
The net diagrams have another remarkable feature, namely that all paths in the diagram that start from 
one nucleus finish at the other nucleus. In other words, there are no loops in the diagram. Each path can 
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be cut at some value of rapidity, and contributes to Eq. (4.6). Keeping only the factor in the exponent we 
obtain that 



Afi (C)/2 = / dzexpl -Z-Z + ) (4.7) 



dz exp 



aji.(y-n/VK-») / a |m(y 7 Vi) 



where z = z + — ^A s V. 

Integrating over 2 using the steepest decent method we obtain that the saddle point in the z-integration 
is z~sp = C/4 which leads to 



Aft (C) /2 = H exp 



JVLU=} ( 4. 



It is interesting to note that the dependence on Y' disappears, and the approach to the unitarity bound 
is much milder than in the case of the solution to the BK equation. Indeed, in the saturation region with 
two saturation scales AQ oc exp ^ — Const ^/Y^j while for the BK equation AO, oc exp ^ — ConstY 2 ^ . 

The function H absorbs all the pre-exponential factors that depend on Y and impact parameters. The 
nucleus-nucleus scattering amplitude can be written using Eq. ( fOID - Eq. (|2~23p in the form 



Iiru4 (Y; B) = 1 - exp ( - J d 2 bd 2 b' (s Al (b) S A2 (b -tj - ) (4.9) 

where b is the impact parameter of the nucleon inside of one nucleus and B is the distance between the 
centers of two nucleons. b' in Eq. ( |4.9| ) is the impact parameter of the dipole inside of the nucleon. As 
we have discussed above, in our approach we integrated over this impact parameter. The typical value of 
b' is our new dimensional parameter, that determines the range of energy which we can reach within this 
approach. Generally speaking the average b' increases with energy, but we assume that < b' > <C Ra- 
This inequality determines the range of energies where we can trust our approach. The energy dependence 
of < b > cannot be found in our approach, since in the framework of perturbative QCD, < b >oc s A . 
However we know (see the Froissart theorem in ref. [35]) that < b' > can increase only logarithmically, but 
the parameters of such an increase will depend crucially on the non-perturbative parameter: the mass of 
the lightest hadron. The best estimate is to insert into Eq. ( fl.9[ ) the experimental value of the inelastic 
cross section for the proton-proton interaction. In doing so we obtain 

W (Y; B) = 1 - exp ( - a in (pp) J d 2 bS Al (b) S A2 [b - ft) ) (4.10) 

It should be stressed that Eq. ( [4. 1Q|) came naturally from our approach, whereas however it is certainly 
incorrect in the usual Glauber- Gribov [36] approximation (see Refs. [37]). 
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5. Conclusions 



In conclusion we can summarize our results for the nucleus-nucleus scattering amplitude as follows. 

1. The contribution of short distances to the opacity 0, dies at high energies. 

2. The opacity tends to unity at high energies. 

3. The main contribution that survives comes from the soft (long distance) processes for large values of 
impact parameters. 



The corrections to the opacity Q = 1 — AN that stem from short distances, have been discussed in this 
paper and it was shown that they behave differently from the corrections to the Balitsky-Kovchegov equa- 
tion. Indeed, it turns out that AN oc exp ^— Const\fY^ while for the BK equation AiV oc exp (— ConstY 2 ^. 

The most salient result of this paper is the formula of Eq. (4.10) that describes nucleus- nucleus colli- 
sions. This formula is instructive, especially since in the usual Glauber-Gribov approach, there is no reason 
to expect that this formula works [37]. It should be mentioned that Eq. ( |4.10|) is used in the interpretation 
of all nucleus-nucleus scattering experiments, including the experiments at RHIC. 

All of these above mentioned results are based on the BFKL Pomeron calculus, which follows from the 
direct sum over Feynman diagrams. However, we should emphasize that the equivalence of this approach 
with other approaches (see refs. [9,12] for example,) that originate from the JIMWLK equation [7], have 
not been proven. It should be stressed that the two strategies reflect two different fundamental features 
of QCD, namely the BFKL Pomeron calculus satisfies, by construction, the t-channel unitarity while the 
JIMWLK-based approaches are precise from the point of view of s-channel unitarity. 
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